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NOTATION 

a,b,c Eq.(64)  coefficients,  defined  by  Eq' s . (65) - (67) 

A laminate  inplane  stiffness  matrix;  Eq.  (34) 

B laminate  coupling  stiffness  matrix;  Eq.  (35) 

C, lamina  stiffness  matrix;  Eq.  (1) 

D laminate  moment  stiffness  matrix;  Eq.  (36) 

e strain 

E  ,E  lamina  longitudinal  and  transverse  stiffnesses 

x*  y  6 

Fvy»Fy>fvv»fv»Fvv»Fq*- • -Tsai-Wu  failure  parameters;  EqTs.(56) 
xx  x   yy   y  xy   ss    and  (5g) 

G   lamina  shear  modulus 

M laminate  moment  vector;  Eq.(26) 

NL laminate  thermal  moment  vector;  Eq.38) 

N laminate  force  vector;  Eq.(25) 

NT laminate  thermal  force  vector;  Eq.(37) 

Q lamina  stiffness  matrix;  Eq.(13) 

R ....lamina  strength  ratio;  Eq.(60) 

S. lamina  compliance  matrix;  Eq.(13) 

t plate  thickness 

T lamina  temperature 

u,v,w displacements  along  the  x,y,z  axes 

x.,y.,z. lamina  axes 

x,y , z. „.... . ...  laminate  axes 

a vector  of  thermal  coefficients  of  expansion 

| vector  of  thermal  coefficients 

Y engineering  shear  strain 

k laminate  curvature  vector 

<{>..., orientation  angle  of  i-th  ply 

9 temperature  of  i-th  ply;  Eq.(2) 

v   ,v   longitudinal  and  transverse  Poisson's  ratio 

xy '   yx  & 

X,X' ,Y,Y' ,S. . . . lamina  strengths 
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1 . 0   Introduction 

In  a  severe  thermal  environment,  materials  generally  ex- 
perience a  degraduation  of  stiffness  and  strength,  and  as  a 
result  there  is  a  diminishment  in  the  load  carrying  capacity 
of  the  material.  In  order  to  design  a  structure  which  will 
perform  safely  in  a  thermal  environment,  the  structural  engi- 
neer must  have  an  analysis  program  which  will  predict  the  be- 
havior of  the  structure.  This  report  documents  an  analysis 
for  the  behavior  of  a  fiber  reinforced  composite  plate  sub- 
jected to  loads  in  a  thermal  environment. 

2. 0   Description  of  Fiber  Reinforced  Composite  Plates 

As  a  result  of  imperfections,  such  as  dislocations  and 
surface  scratches,  the  strength  of  bulk  material  is  orders  of 
magnitude  less  than  the  strength  of  the  material  in  crystal 
or  whisker  form.  The  cost  of  producing  crystals  is  so  great 
as  to  prohibit  production  except  on  a  limited  basis.  On  the 
other  hand,  fibers  of  a  material  can  be  produced  at  a  cost 
which  permits  their  application  to  structural  systems. 

Depending  upon  it's  diameter,  a  fiber  of  a  material  has 
a  strength  intermediate  to  the  bulk  and  crystal  strengths  of 
the  material.  For  example,  the  strength  of  steel  fiber  is 
approximately  4.14  MPa  (600,000  psi) ,  compared  to  a  strength 
of  1.04  MPa  (150,000  psi)  for  steel  in  it's  common  bulk  form. 
Thus  the  use  of  steel  fiber  in  place  of  common  steel  results 
in  a  four-fold  reduction  in  weight.  Actually  the  use  of  fiber 
requires  their  embedment   in  a  matrix  of  another  material  in 
order  to  achieve  an  useful  structural  form.  However  even  with 
embedment,  fiber  reinforced  materials  offer  a  significant  in- 
crease in  the  strength  to  weight  ratio  over  common  materials. 
This  is  the  primary  reason  for  using  composite  materials. 

One  of  the  most  common  useful  structural  component  is 
the  plate  element.  A  fiber  reinforced  composite  plate  (or 
simply,  a  laminate)  consists  of  layers  of  lamina  (or  plies) 


stacked  one  upon  another.  A  lamina  itself  is  comprised  of  uni- 
directional fibers  embedded  in  a  binding  material,  called  the 
matrix.  Generally,  the  binding  matrix  material  is  a  much  weak- 
er and  more  ductile  material  than  the  fiber  reinforcement  em- 
bedded within  it.  Some  common  fiber/matrix  combinations  are 
graphite/epoxy ,  glass/epoxy,  and  boron/aluminum.  Some  idea  of 
the  dimensions  involved  can  be  obtained  from  the  dimensions 
of  a  typical  graphite/epoxy  lamina; 

i)  fiber  diameter  is  approximately  .008  mm.  (.0003  in.) 
ii).  lamina  thickness  is  approximately  .125  mm.  (.005  in.) 
iii)  fiber  volume  to  total  volume  (called  the  fiber 
fraction)  is  about  0.5 
For  these  dimensions,  a  12.5  mm.  (0.5  in.)  thick  laminate 
will  consist  of  100  plies  (lamina)  stacked  one  upon  another. 

As  a  result  of  unidirectional  fibers  embedded  within  it, 
a  lamina  is  not  isotropic.  The  stiffness  of  a  lamina  along  it's 
fiber  direction  (called  the  longitudinal  stiffness)  is  much 
greater  than  the  stiffness  of  the  lamina  in  the  direction 
orthogonal  to  the  fiber  (called  the  transverse  stiffness) . 
For  example,  in  the  case  of  a  typical  graphite/epoxy  lamina, 
the  longitudinal  stiffness  (E   or  E^)  and  the  transverse 

stiffness  (E  or  ET)  are  181  GPa  (26xl06  psi)  and  10.3  GPa 

6    ' 

(1.5x10   psi),  respectively. 

Stacking  lamina  of  different  orientation   and  thick- 
ness upon  one  another  to  form  a  composite  plate  results  in  an 
orthotropic  body.  One  of  the  most  attractive  features  of  a 
composite  plate  is  the  capacity  for  design  of  material  prop- 
erties, an  option  not  previously  available  to  structural 
engineers . 

3. 0   Lamina  and  Laminate  Coordinate  Systems 

In  dealing  with  fiber  reinforced  composite  plates,  it  is 
necessary  to  employ  two  types  of  coordinate  systems.  The  first 
type  of  coordinate  system  is  associated  with  a  lamina.  This 
coordinate  system  is  referred  to  as  the  natural  coordinate 


system,  or  the  local  coordinate  system,  or  the  lamina  (or  ply) 
coordinate  system.  Here  we  assign  the  cartesian  pair  of  x,y 
axes  to  the  plane  of  the  lamina,  with  the  x-axis  along  the 
fiber  direction  of  the  lamina.  Since  adjacent  lamina  have 
different  fiber  orientation  with  respect  to  one  another,  there 
are  as  many  natural  coordinate  systems  as  there  are  plies.  The 
i-th  lamina  (ply)  has  the  coordinate  system  (x,y) .  or  (x. ,y.) . 
This  is  shown  in  Fig.l  below. 
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Figure  1.  Lamina  Layup  to  Form  a  Laminate 


The  second  type  of  coordinate  system  associated  with  a 
composite  plate  is  the  coordinate  system  which  ties  all  of 
the  lamina  coordinate  systems  to  a  single  pair  of  axes, 
called  the  laminate  or  global  coordinate  system,  and  denoted 
by  (x,y) .  These  axes  lie  in  a  plane  parallel  to  the  lamina 
planes,  and  are  generally  located  at  the  midplane  of  the  la- 
minate, as  shown  in  Fig. 2  for  a  3  ply  laminate. 

Each  of  the  local  (lamina)  coordinate  systems  are  rela- 
ted to  the  global  (laminate)  coordinate  system  by  the  angle 
<J>.,  where  $.  is  the  angle  between  the  x.  axis  of  the  i-th 
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lamina  and  the  x  axis  of  the  laminate.  In  this  way,  the 
orientation  layup  of  a  n-ply  laminate  can  be  denoted  by  the 
notation  [<t>  1  ,  $2  ,  4>3  ,....,  <j>]  .  There  is  no  requirement  that 
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Figure  2.  Lamina  and  Laminate  Coordinate  Systems 

the  lamina  be  of  equal  thickness,  although  this  is  a  common 
construction. 


4. 0   Stress-Strain  Relations  in  Natural  Coordinates  for  a 

Lamina  in  Plane  Stress 

In  the  case  of  plates  subjected  to  inplane  and  trans- 
verse loads,  the  result  is  the  plane  stress  state.  This  means 
that  for  a  lamina  in  the  x,y  plane,  interlaminar  stresses  a   , 

zy'      zx  are  negligible  compared  to  inplane  stresses  cr   , 

a   ,  and  a   ,  and  are  therefore  ignored.  It  should  be  noted 

yy '      xy '  6 

in  passing  that  interlaminar  stresses  can  be  very  large  in 
local  regions  of  boundaries.  (See  Ref.l  for  details).  In  this 
report  interlaminar  stresses  are  not  considered. 
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In  terms  of  it's  own  natural  coordinate  system,  the 
stress-strain  relations  for  the  i-th  lamina,  in  tensor  nota- 
tion, are 

(a.,).  =  (C,   e   -  3-i  6)-    j ,k,r,s=x- ,y.      (1) 
v  jk'i  *■  jkrs  rs   jk  Ji        J '    *    '        i  'i 

where  a.-,    and  e   are  the  second  order  stress  and  strain 
tensors,  3-^  is  the  second  order  tensor  of  thermal  coeffi- 
cients, C.,    is  the  fourth  order  tensor  of  stiffness 
j  krs 

coefficients,  and  0  is  the  difference  in  temperature  from 
the  stress  free  temperature.  That  is,  if  T  is  the  stress  free 
(curing)  temperature,  and  T.  is  the  temperature  of  the  i-th 
lamina,  then 

ei-(Ti-  T)  (2) 

Throughout  this  report,  the  Einstein  summation  convention 
is  in  effect  for  all  tensor  equations. 

Alternatively,  Eq.(l)  may  be  written  for  strains  in 
terms  of  stress,  as 

(ejkV(Sjkrs°rs+  ajk9)i   J.k.r,.-*^     (3) 
where  S.,    is  the  fourth  order  tensor  of  compliance  coeff- 
icients, and  a. k  is  the  tensor  of  thermal  expansion 
coefficients.  We  note  that  properties  of  a  lamina,  such  as 
stiffnesses,  compliances,  thermal  coefficients,  and  thermal 
expansion  coefficients,  are  functions  of  temperature.  How- 
ever over  narrow  ranges  of  temperature,  properties  can  be 
treated  as  constant  without  introducing  significant  error 
into  the  analysis. 

As  a  matter  of  convenience,  hereafter  we  proceed 
using  contracted  "vector"  notation.  It  should  always  be 
kept  in  mind  that  stresses,  strains,   thermal  coefficients, 
thermal  expansion  coefficients,  stiffnesses,  and  compliances 
are,  in  fact,  tensors,  not  vectors,  and  therefore  obey  the 
appropriate  tensor  transformation  laws  when  transformations 
occur.  The  second  order  tensors  (a .  . ,e-  .  ,  3.. ,  and  «..)  in 
the  contracted  notation  have  a  single  subscript  in  accord- 
ance with  Table  1,  and  the  fourth  order  tensors  (S- -^i  anc* 
C.,,)  carry  the  double  subscript  notation  shown  in  Table  2. 


Tensor 

Contracted 

Subscripts 

Notation 

xxxx 

11 

xxyy 

xxxy  and  xxyx 

12 

Tensor 

Contracted 

16 
21 
22 

Subscripts 

Notation 

yyxx 

yyyy 

yyxy  and  yyyx 

XX 

1 

26 

yy 

xy  and  yx 

2 

xyxx  and  yxxx 

61 

6 

xyyy  and  yxyy 
xyxy  and  yxyx 

62 
66 

Table  1. Contracted  Notation   Table  2. Contracted  Notation  for 
for  2   order  Tensors         4   order  Tensors 


For  example,  using  the  symmetry  conditions,  e   =e   , 

C-v-yy-w=  Cyyvv  an<*  s0  on»  anc*  defining  the  engineering  shear 

strain  Y   =(e   +  e   ),  equations  (1)  become 
Ax   Ay   yA 

[a1  =  Cnei  +  C12e2  +  C16y6  -  3^ 

[a2  =  C2ie;L  +  C22e2  +  C26y6  -  020]i 

ta6  "  C61el  +  C62e2  +  C66^6  '  ^i 


(4) 


or,  in  vector  notation 

ai=  (  Ce  -  36  )i  (5) 

where  the  i  subscript  denotes  the  equation  is  for  the  i-th 
lamina.  In  equation  (5),  and  hereafter,  a  single  tilda  under 
a  quantity  denotes  the  quantity  is  treated  as  a  3x1  vector, 
a  double  tilda  under  a  quantity  means  the  quantity  is  treated 
as  a  3x3  matrix,  and  a  quantity  without  a  tilda  beneath  it 
denotes  that  the  quantity  is  a  scalar. 

As  a  result  of  the  natural  coordinates  being  along  and 
orthogonal  to  the  fibers  of  a  lamina,  the  natural  coordinates 
are  the  principal  axes  of  the  lamina.  Thus  the  normal  and 
shear  effects  uncouple,  and  Eq.  (5)  takes  on  the  special  form, 


I 


i   - 


t . 

1 


11 


21 


'12 


'22 


66 


f      s 

el 


c2 
Y6 


,  3. 


(6) 


That  is,  in  the  natural  coordinates  of  a  lamina,  C,.=C0£=C, 

'  '   16   26   66 

=0,  and  36=0.  Furthermore,  as  the  xy  plane  is  a  plane  of 
symmetry,  ^i^'^ll*    anc*  therefore  C  is  a  symmetric  matrix. 
Thus  a  lamina  has  four  independent  stiffness  coefficients 
(C-,  -,  ,C,  2  »C72 »  anc*-  ^66^  '  an(*  two  independent  thermal  coeffi- 
cients (3-t  ,  and  3?)  •  Such  materials  are  called  orthotropic 
materials.  It  is  important  to  note  that  the  special  form  of 
C  in  Eq.(6),  associated  with  uncoupled  normal  and  shear  be- 
havior, is  valid  only  in  natural  coordinates  of  a  lamina. 

The  notation  used  most  frequently  in  the  composites 
literature  is  §  for  compliance,  and  0  rather  than  C  for 
stiffness.  Hereafter  we  shall  use  S  and  Q.  In  contracted 
notation,  Eq's.(l)  and  (3)  are, 


and 


2i  =  (?!  "£e)i 
e.  =  (So  +a6). 


(7) 
(8) 


Premultiplication  of  Eq.(7)  by  Q"   and  solving  for  strain 
e  gives 


Si  =  eg'1?   +  Q"1^) 


(9) 


A  comparison  of  equations  (8)  and  (9)  shows  the  relations 
between  Q  and  S,  and  3  and  a  are 

,-1 


and 


2  =  i 

9  =  Q"1?  • 


(10) 

(ii) 


The  stiffness  coefficients  Q . ,  are  related  to  the  common 

3  K 


engineering  coefficients  E  ,E v,vvv,  and  G   ,  by  Eq's.(12) 
The  engineering  coefficients  are  obtained  by  performing  the 
simple  tension  and  shear  tests  on  a  lamina  as  depicted  in 


Figure  3.  Pulling  a  lamina  along  the  fiber  direction  (Fig. 3a) 

gives  the  longitudinal  stiffness,  E  .  The  tension  test  in  the 

orthogonal  direction  (Fig. 3b)  gives  the  transverse  stiffness 

E  .  The  longitudinal  Poisson's  ratio  v    is  the  ratio  of  e 
y  xy  x 

to  e   obtained  in  the  longitudinal  tension  test.  The  long- 
itudinal shear  modulus  G   is  obtained  from  the  shear  test 

xy 

shown  in  Fig. 3c.  The  relations  between  S.,  and  engineering 
coefficients  are, 


Sn  =  (l/Ex) 
S22=(l/Ey) 


(12) 


Note  that  the  second  of  these  relations  defines  a  transverse 
Poisson's  ratio,  v  .  The  x  and  y  subscripts  denoting  longi- 
tudinal and  transverse  properties  are  quite  often  replaced  by 
T  and  L  subscripts,  i.e.,  EL  denotes  the  longitudinal  stiff- 
ness of  a  lamina,  ET  is  the  transverse  stiffness,  and  so  on. 
Using  Eq.  (10)  ,  the  relations  between  stiffness  coefficients 
Q . •  and  the  engineering  constants  are  found  to  be, 

Ql2=(521=CvxyEx/t1-vxyV))  =  CW(1"VxyV:i) 
Q22=(V(1-Vvx)  (13) 


Q66~Gxy 


Figure  3.  Experimental  Tests  for  Engineering  Stiffness  Constants 
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5. 0   Stress-Strain  Relations  for  a  Lamina  in  Plane  Stress 

with  Respect  to  Laminate  Axes 

In  order  to  develope  the  equations  for  a  laminate,  it  is 
necessary  to  relate  each  lamina  comprising  the  laminate  to  a 
common  coordinate  system,  the  laminate  axes.  This  is  accomp- 
lished by  coordinate  transformation  from  lamina  axes  to  lami- 
nate axes.  Each  lamina  coordinate  system  (x,y,z).  is  oriented 
to  the  laminate  axes  (x,y,z),  where  z.  and  z  coincide,  by  an 
angle  <j>..  Recall  from  Fig.  2,  that  <p .  is  the  angle  between  the 
x.  axis  of  the  i-th  lamina  and  the  x  axis  of  the  laminate. 

Hereafter  all  quantities  with  a  bar  over  the  quantity 
symbol  means  that  the  quantity  is  with  respect  to  laminate 
axes,  and  quantities  without  a  bar  over  them  means  the  quant- 
ities are  with  respect  to  lamina  axes.  Thus  the  stresses, 
strains,  thermal  coefficients,  and  stiffnesses  with  respect 
to  the  laminate  coordinate  system  are  a,e,£,  and  Q.  Using  the 
tensor  transformation  laws  for  the  second  order  tensors, 

(ajk}i  =  (1jm1knamn)i 

(e.  J .  =  (1.  1,  e   ) .  (14) 

CBjk}i  =  (1jm1knUi 

where  1    is  the  cosine  of  the  angle  between  the  r-th  lamina 
axis  and  the  s-th  laminate  axis,  gives  the  expressions  for 
lamina  quantities  with  respect  to  laminate  axes.  Denoting 
the  transformation  between  coordinate  systems  by  the  matrix 
T,  Eq's.(14)  may  be  written, 

o.  =  (Ta) . 

§i  =  (Ii}i  (is) 


Substitution  of  Eq's.(15)  into  Eq's.(7)  gives 

T . a .  =  Q.T.e.  -  T  .  3  .  6  . 
~i~i   £i*i~i   *i~i  1 

Premultiplication  through  by  the  inverse  of  T.  gives, 
a.    =  TT1Q.T.e.  -  8.6. 

Denoting  the  matrix  T.  Q.T.  by  Q.,  where  Q.  is  the  3x3  stiff - 

« l  ~i«i    ~i        3i 

ness  matrix  of  the  i-th  lamina  with  respect  to  laminate  axes, 
we  obtain  the  stress-strain  relations  for  a  lamina  in  terms 
of  laminate  axes  as, 

a •=  0-e   -  §  -  8  . 
~i   5i5i    £i  i 

Recalling  Eq.(ll),  yields  the  more  useful  expression  (since 
the  coefficients  of  thermal  expansion  are  generally  known) , 


a.=  Q. (e.-  o.e.) 


(16) 


where  the  i  subscript  refers  to  the  i-th  lamina.  Explicitely 
Eq' s . (16)  are, 


LC6J 


;11 


12 


:16 


§21   §22   §26 
§61   §62   §66 


'   ~eI 

\ 

~e2 

- 

e. 

i 

a2 

i 

L^6. 

i 

ft. 

(17) 


where,    upon  using   c.=   cos<J>.    and   s-=    sin^.,    we  obtain, 

(Qn)i=  [Qnc4   ♦    2(Q12  +  2Q66)s2c2    ♦   Q„.4]i 

(Q12)i=  [CQu*Q22-4Q66)s2c2   ♦   Qu(s4+c4)]i 

(Q22)i=  [Qn=4  +2(Q12+2Q66;)s2c2    ♦  q22c4]i 

CQieJ±-  [CQn-Q12-2Q66)sc3   ♦    (Q12-Q22  +  2Q66)s3c]1 

(Q26)i=  [CQn-Q12-2Q66)s3c   ♦    CQ12-Q22*2Q66^sc3]i 

®66h~  [CQu+Q22-2Qi2-2Q66^s2c2  +  Q66^4+c4^i 


(18) 
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We  note  from  Eq's.(18)  that  the  stiffness  matrix  of  a 
lamina  with  respect  to  the  laminate  coordinate  system  is  a 
full  3x3  matrix.  However,  there  are  still  only  four  inde- 
pendent coefficients;  any  two  of  the  Q-,  coefficients  may 

J  K 
be  written  as  functions  of  the  remaining  four  coefficients. 

The  Q  matrix  is  symmetric,  that  is,  Q,2  =  Q?i»  §16  =  §61* 

and  §26  =  Q62- 

The  fact  that  Q  is, in  general,  full,  means  that  coupling 

between  normal  and  shear  behavior  occurs.  In  contrast,  for 
isotropic  materials,  normal  and  shear  behaviors  uncouple. 
For  cross  ply  laminates,  that  is  laminates  with  orthogonal 
plies,  subjected  to  inplane  forces  and  bending  moments, 
normal  and  shear  behaviors  uncouple.  For  angle  ply  laminates, 
that  is  laminates  with  some  non-orthogonal  plies,  coupling 
between  normal  and  shear  behaviors  occurs.  This  difference 
in  behavior  between  cross  ply  laminates  and  angle  ply  lami- 
nates is  shown  by  an  illustrative  example  in  section  11. 

6 . 0   Strain-Displacement  Relations  for  a  Laminate 

In  accordance  with  common  thin  plate  theory,  the  as- 
sumption that  planes  before  loading  remain  planes  after  load- 
ing (Kirchhoff -Love  theory)  is  taken  as  the  basis  of  deform- 
ation of  a  laminate.  In  this  case,  the  displacement  field  is 
given  by, 

ii  =  u   -  z(3w/3x) 

(19) 


v  =  vQ  -  z(3w/3y) 


} 


where  u  and  v  are  the  displacements  in  the  x  and  y  directions, 

ii  and  v   are  the  midplane  displacements,  w  is  the  displace- 

o      o   _  _ 

ment  in  the  z  direction,  and  z  is  the  distance  from  the  mid- 
plane  to  the  point  of  interest. 

Using  the  linear  strain-displacement  relations, 


} 


i   =  (3u/3x)      e   =  (3v/3y) 
Yvv  =  Ou/3y)  +  (3v/3x) 
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equations  (19)  become, 

e±   =  (3u0/3x)-  z^w/Sx2) 

e2  =  (3vQ/3y)-  z^w/Sy2)  ^     (21) 

Y6  =  (3uQ/3y  +  3v  /3x)  -  2z  (^w/SxSy) 

Equations  (21)  may  be  written  in  vector  form  as, 

e  =  eQ  +  zk  (22) 

where  e   is  the  strain  vector  associated  with  midplane  exten- 
sional  displacements  given  by 

Ce1)Q  3uQ/3x  S 

(e2)Q   =   3vo/3y  \  (23) 

^6}o      8Go/9^  +  3V3*   / 

and  k  is  the  curvature  vector  associated  with  bending  of  the 
midplane,  given  by 

(kx)     -  ?vfizx  \ 

(k2)  =   -  tfutfy  \  (24) 

(k6)    -232w/3x3y       ) 

Since  the  midplane  strains  and  curvatures  belong  to  the  lami- 
nate, they  do  not  change  from  one  lamina  to  another,  and 
therefore  the  bar  over  these  quantities  is  not  required  and 
henceforth  will  be  dispensed  with. 

7 . 0   Force-Displacement  Relations  for  a  Laminate 

Equation  (22)  shows  that  the  strain  at  a  point  is  the 
sum  of  two  components,  (i)  midplane  extensional  strain,  and 
(ii)  midplane  curvature.  Figure  4.  depicts  the  displacement, 
strain,  modulus,  and  stress  field  with  respect  to  the  x,z 
plane  at  an  arbitrary  y  location.  The  displacement  and 
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strain  fields  are  shown  together  in  Fig. 4b,  while  in  fact, 
one  is  a  scalar  times  the  other.   In  accordance  with  the 
Kirchhof f -Love  assumption,  the  displacement  field  is  linear 
across  the  laminate  thickness.  As  a  result  of  the  jump  dis- 
continuity of  modulus  from  lamina  to  lamina  (Fig. 4c),  the 
stress  field  also  has  a  jump  discontinuity  from  lamina  to 
lamina  (Fig.4d) . 


X 


Laminate 
Field 

(a) 


Displacement  (and 
Strain)  Field 

(b) 


Modulus  Field   Stress  Field 


(c) 


(d) 


Figure  4.  Displacement,  Strain,  Modulus  and  Stress  Distribu- 
tions 


The  laminate  forces  N  and  moments  M  associated  with  a  la- 
minate are  obtained  by  integrations  of  the  stress  field  over 
the  laminate  thickness.  The  extensional   laminate  forces  are 
given  by, 


r  -1 
fi] 

N.. 

N, 


f 


-t 
2 


dz 


or   N   =/  adz 


■/ 


(25) 
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and  the  laminate  moments  are  given  by, 


zdz 


1 

ftl 

-t 

r 

u.    6- 

8-/2 


zdz    (26) 


where  t  is  the  plate  thickness.  The  laminate  forces  and  mo- 
ments are  shown  in  Figure  5.  Components  N,  and  N2  are  the 
normal  components  of  force,  and  N6  is  the  shear  component.  M, 

and  M2  are  bending  moments,  and  M^  is  the  twisting  moment. 


*.y,** 


Figure  5.  Laminate  Forces  and  Moments 


Since  the  stress  field  across  the  laminate  thickness  is 

the  sum  of  the  lamina  stresses, 

n 

a  =  /   a . 

-  fcf1 

equations  (25)  and  (26)  become 

.z. 


and 


N  =>   /  a. 
3.-1 

-      #    f-H 

M  =  >  /  a.  z< 


dz 


(27) 


(28) 


zi-l 
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where  z . _ ,  and  z.  are  the  upper  and  lower  limits  of  integra- 
tion for  the  i-th  lamina,  and  n  is  the  number  of  lamina  com- 
prising the  laminate. 

Substituting  the  expressions  for  lamina  stresses  in 
terms  of  strains,  Eqrs.(16),  into  the  previous  Eq's.(27)  and 
(28)  gives  the  laminate  forces  and  moments  in  terms  of  strains 
as , 

N  =   §  j    ^iCgi"  <*iei)dz  (29) 

i=1z-  ? 
and  Zi-1 

8  -  ?  /  iQiCe±-  ai8i)£d£  (30) 

i-lii-f 

The  strain  in  the  i-th  lamina  is  given  by  Eq.(22), 

e.  ■  e  +  z  k  (31) 

where  e  and  k  are  the  midplane  extensional  strain  and  mid- 
~o  r 

plane  curvature  vectors  respectively,  and  are  constant  from 
one  lamina  to  another.  Substituting  Eq.(31)  into  Eq's.(29) 

and  (30)  gives, 

z . 

dz 


N  =  I  f0.(e.  +  zk  -  a. 9.) 
i=rz.  *x  °         1   1 


M  =  t    f  Q-(e   +  zk  -  a.9.)zdz 

•_J   il  -O      ~  ~1  1; 

i-lzi."! 

Since  e   and  k  do  not  vary  with  z,  and  Q.,a. ,  and  9.  are  con- 

~o                  '  Si  ~1        1 

stant  for  a  lamina,  the  previous  equations  become, 

n             Zi  z 


N  =r  Q-[(e  -a. 9.)  fdz  +  ic  f  zdz] 
.  ,  Si  ~o  ~i  i  J      ~  J 


i-1    zi-l 


n        .     zi.  _  .  zi_  _ 

M  =  Z  Qi[(eQ-ai9i)  /  zdz  +  ic  J  zdz] 

i=l  s    "*        z .  ,  z .  , 

l-l  l-l 


Upon  integration,  we  obtain, 


N  =  Z     QiCe0-ai9i)(£i-z._1)  +  ?  Q^^2!' zi-l)  (32) 
i=l  ~  i=l  z  2 
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and 


Z  2  Z  2    -,      „  3    3 


ft-   Z   Q.(e  -a.e.) (zi"zi-l^  +  z  Q.k.(z.-z.  J 
i=l  z1   ~°   -1  x  — i=l  -J-^J-   J-   i-l 


Defining  the  matrices, 

n 


A  =  z  Qi(Vzi-l}  C34) 

~  i»l  = 

B  =  z     Q.  (z?-Z?  J/2  (35) 

1  =  1 

D  =  ?  Q^q-z^^/e  (36) 

s  i=l  = 


and  vectors, 


n 


5i  \h  9i2iWxi-i'  <"> 

h  Vi  3i2i9i("zr£!-i)/2  .  (38) 

equations  (32)  and  (33)  become, 

N  =  Ae   +  Bk  -  NT  (39) 

M  =  Be   +  Dk  -  MT  (40) 

Equations  (39)  and  (40)  are  the  "force-displacement"  rela- 
tions for  a  laminate. 

The  NT  and  MT  vectors  are  the  so  called  thermal  force 
vector  and  thermal  moment  vector  respectively,  as  they  result 
from  a  thermal  environment.  If  the  temperature  of  each  lami- 
na is  equal  to  the  stress  free  temperature  (i.e.,  the  curing 
temperature),  then  9.  is  zero  for  each  lamina,  and  the  therm- 
al vectors  are  equal  to  the  zero  vector.  In  this  case,  Eq's. 
(39)  and  (40)  reduce  to  the  force-displacement  relations  for 
a  laminate  not  subjected  to  thermal  effects. 

The  A  matrix  relates  inplane  forces  to  inplane  strains, 
and  the  D  matrix  relates  moments  to  curvature.  The  B  matrix 
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couples  inplane  forces  to  curvature,  as  well  as  moments  to 
inplane  strains.  In  the  special,  albeit  common,  case  of  sym- 
metric laminates,  the  contribution  to  S  from  a  lamina  above 
the  midplane  is  cancelled  by  the  contribution  from  a  counter- 
part lamina  below  the  midplane,  and  as  a  result,  the  B  matrix 
is  the  zero  matrix,  i.e.,  B  =  0.  Thus  for  symmetric  laminates, 

ST  ■  £§0-  5t  (4i) 

MT  =  Dk  -  MT  (42) 

and  therefore   the  extensional  and  bending  behaviors  are  un- 
coupled. If  the  laminate  is  symmetric  and  the  temperature 
distribution  is  symmetric  with  respect  to  the  midplane,  then 
from  Eq.(38),  the  thermal  moment  vector  is  equal  to  the  zero 
vector  and  Eq. (42)  reduces  to 

M  =  Dk  (43) 


8 . 0  Stress-Force  Relations  for  a  Lamina  in  a  Symmetric 
Laminate 
In  order   to  determine  whether  a  ply  has  failed,  it  is 
necessary  to  determine  the  stresses  in  each  ply.  The  devel- 
opement  of  the  stress-force  relations  for  a  lamina,  which 
follows,  is  restricted  to  symmetric  laminates.  This  is  not  a 
severe  restriction  since  symmetric  laminates  are  in  most 
common  use.  Solving  Eq's.(41)  and  (42)  for  strains  and  curv- 
atures gives, 

e   =  A_1(N  +  NT)  (44) 

and  -, 

K   =  D_i(M  +  M~)  (45) 

~     £    ~    ~  i 

The  lamina  strains  are  related  to  the  midplane  strains 
and  curvatures  of  the  laminate  by  Eq. (31) , 

e.  =  e  +  zk  (31) 

-i   ~o 
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Substitution  of  Eq's.(44)  and  (45)  into  Eq.(31)  gives, 

ei  =  A_1(N  +  NT)  +  ijD"1^  +  Kj.)  (46) 

where  z.  is  the  distance  from  the  laminate  midplane  to  the 
middle  of  the  i-th  lamina,  that  is, 

z.  =  (z.  +  z.  J/2  (47) 

l     l    l-l/'  v  } 

Lamina  stresses  with  respect  to  laminate  axes  are  related  to 
lamina  strains  in  laminate  axes  by  Eq.(16)  as, 

a.  =  Q. (e.  -  i.e.)  (16) 

~i   ~i  ~i   ~ii  K      J 

Finally  then,  substitution  of  Eq.(46)  into  Eq.(16)  gives  the 
lamina  stresses,  with  respect  to  laminate  axes,  in  terms  of 
laminate  inplane  forces  N,  laminate  moments  M,  thermal  forces 
NL ,  and  thermal  moments  MT  as, 

ai=  QitA'^N  +  NT)  +  ZiP"1^  +  ftj.)  -§i6i]  (48) 


Equation  (48)  takes  on  special  forms  for  particular  cases;  a 
few  of  which  are  presented  below. 

a)  Temperature  free  effects  occur  when  the  plate  is  at  a 
uniform  temperature  equal  to  the  curing  temperature  of 
laminate.  In  this  case,  9.  =0,  and  therefore  NT  and  MT 
are  equal  to  the  zero  vector,  and  Eq.(48)  reduces  to, 

a.  =  Q. (A_1N  +  z.D_1M)  (49) 

In  the  absence  of  lateral  loads  on  the  laminate,  M  =  0, 
and  Eq.(49)  becomes, 

a   =  QA_1N  (50) 

b)  Absence  of  moments  occurs  when  the  laminate  is  not  sub- 
jected to  lateral  loads.  In  this  case,  M  =  0,  and 
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equation  (48)  reduces  to, 

i±    =  SitA^CN  +  NT)  +  ZiP"1^  -ai9i]        (51) 

If  the  temperature  field  is  symmetric  with  respect  to 
the  midplane  of  the  laminate,  then  $L,   is  the  zero 
vector  and  Eq.(51)  becomes, 

a.  =  Q.  [A_:L(N  +  NT)  -  a. 9.]  (52) 

c)  Thermal  effects  only  occur  when  the  laminate  is  not 
subjected  to  inplane  and  lateral  loads.  In  this  case 
N  =  M  =  0,  and  Eq.(48)  reduces  to 

2i  =  §i[  h'\   +  ^iS"1?!  "  ^iei]        (53) 

If  the  temperature  field  is  symmetric  with  respect  to 
the  laminate  midplane,  then  M^  =  0,  and  Eq.(53)  gives, 

a.    =  Q. [A_1NT  -  a.9.]  (54) 


9 . 0   Failure  Criterion 

When  a  composite  plate  is  subjected  to  a  load  and/or 
thermal  environment,  stresses  develope  in  each  lamina  in  ac- 
cordance with  Eq.(48).  Depending  upon  the  magnitude  of  the 
lamina  stresses,  one  or  more  lamina  may  fail.  When  a  lamina 
fails,  there  are  two  effects.  First,  if  we  assume  brittle 
failure,  then  a  lamina  will  relinquish  it's  load  carrying 
capacity  upon  failure,  and  there  will  be  a  redistribution 
of  load  to  the  remaining  non  failed  lamina.  Secondly,  when  a 
lamina  fails,  the  stiffness  matrices  A  and  D  change   in  ac- 
cordance with  Eqfs.(34)  and  (36),  and  thus  the  redistribu- 
tion of  stresses  in  the  non  failed  lamina  is  not  simply  a 
change  in  magnitude,  but  rather  a  change  in  the  stress  state. 

The  hypothesis  of  a  brittle  failure  of  a  lamina  is  an 


19 


idealization  of  the  actual  behavior.  Actually,  lamina  have 
some  ductility  and  therefore  continue  to  carry  some  fraction 
of  load,  and  contribute  some  stiffness  to  the  laminate  upon 
failure.  Here  we  assume  a  failed  lamina  serves  no  useful  pur- 
pose after  it's  failure. 

As  will  be  demonstrated  by  illustrative  examples  in 
section  11,  failure  of  a  ply  does  not  necessarily  mean  that 
the  laminate  has  failed.  There  we  will  show  that  a  cross  ply 
laminate  subjected  to  uniaxial  loading  does  not  fail  upon 
failure  of  the  first  ply,  and  therefore  the  magnitude  of  the 
load  can  be  increased  beyond  first  ply  failure.  We  also  show 
that  under  uniaxial  loading,  angle  ply  laminates  fail  when 
the  first  ply  fails. 

There  is  no  way  to  establish  the  validity  of  any  failure 
criterion  on  theoretical  grounds.  All  that  can  be  done  is  to 
show  that  a  failure  criterion  gives  analytical  results  which 
agree,  for  a  wide  class  of  problems,  with  experimental  re- 
sults. There  are  several  failure  theories  which  give  good  re- 
sults and  are  in  common  use.  The  most  general  of  the  macro- 
scopic failure  criteria,  which  includes  many  of  the  other 
criteria  as  special  cases,  is  the  Tsai-Wu  criterion,  Ref.(2). 
It  is  the  failure  criterion  adopted  in  this  work. 

According  to  the  Tsai-Wu  criterion,  a  lamina  fails  when 
the  general  quadratic  equation  (in  tensor  form) , 

F-.a-a.  +  F.a.  =  1     i,j=l,2,,,6     (55) 
ij  l  ]     11  >j»>>> 

is  satisfied.  In  Eq.(55),  F.  and  F-.  are  strength  tensors  of 
second  and  fourth  order,  and  the  equation  is  in  contracted 
notation.  In  the  case  of  plane  stress,  the  stress  state  for 
composite  plates,  Eq.(55)  in  terms  of  natural  coordinates  of 
a  lamina,  becomes 

Fllai  +  2F12ala2   +  F22a2  +  F66a6  +  Flal  +  F2a2  =  X  <56) 
Equation  (56)  is  applied  to  each  lamina.  If  the  left  hand 
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side  of  Eq. (56)  is  less  than  unity,  then  the  lamina  has  not 
failed.  If  the  left  hand  side  is  greater  than  unity,  then  the 
lamina  has  failed.  Failure  of  a  lamina  is  imminent  (impend- 
ing) when  the  left  hand  side  is  equal  to  unity,  i.e.,  when 
Eq.(56)  is  satisfied  identically. 

In  order  to  use  Eq.(56)  directly,  that  is  without  modi- 
fication, it  is  necessary  to  use  lamina  stresses  with  respect 
to  lamina  axes.  The  use  of  lamina  stresses  with  respect  to 
laminate  axes  would  require  the  transformation  of  the  second 
and  fourth  order  tensors,  F.  and  F-.,  to  laminate  axes.  Al- 
though this  could  be  done,  it  is  more  expeditious  to  trans- 
form the  lamina  stresses  with  respect  to  laminate  axes, 
given  by  Eq's.(48),  into  lamina  stresses  with  respect  to  la- 
mina axes.  The  transformation  is  accomplished  by  the  use  of 
Eq.(14).  Thus,  we  obtain, 

ia1)i   =  C^.cf  ♦  (52)iS2  ♦  2Co6)icisi 

Ca2).  =  C51)is|  ♦  (52).c|  -  2(56)icisi       (57) 

(c^)  •  ■■(ff1)cs.  +  (o-)  -c  .s  .  +  (ac)  •  (c?-s?) 
*-  6' 1    v  1J    l  i  *•  2'i  i  i  v  6'iv  l  x 

where  c.  =  cos<j>.,  s.  =  sin<{>.,  and  i  denotes  the  lamina. 
i       i*   i      ri' 

Five  of  the  six  strength  parameters  of  Eq. (56)  may  be 
determined  by  a  series  of  simple  experiments  as  follows.  A 
tension  test  of  a  lamina  along  it's  fiber  direction  gives 
the  longitudinal  tension  strength,  X.  A  compression  test 
along  the  fiber  axis  gives  the  longitudinal  compressive 
strength,  X' .  Tension  and  compressive  tests  in  the  transverse 
direction  (i.e.,  orthogonal  to  the  fiber),  gives  the  trans- 
verse tension  and  compressive  strengths,  Y  and  Y' ,  respect- 
ively. A  shear  test  of  the  lamina  gives  it's  shear  strength, 
S.  The  strength  parameters  in  Eq.(56)  are  related  to  the 
experimentally  determined  strengths  as  follows.  Applying 
Eq. (56)  to  each  of  the  five  stress  states  associated  with 
the  five  tests  gives  five  simultaneous  equations  which  may 
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be  solved  for  F,  ,  >F22  ,F66  ,F1 '    an<*  F2   For  examPle»  Eq.(56) 
for  the  longitudinal  tension  test  becomes, 

F11X2  +  F1X  =  -1 
and  for  the  longitudinal  compression  test,  Eq.(56)  is, 

F11X'2"  Fix'=  l 
Solving  these  equations  for  F,,  and  F,  gives  the  first  two 
of  Eq's.(58).  The  remaining  strength  parameters  are  obtained 
in  a  similar  way. 

Fu  =  Cl/XX') 

F2   =  (1/X)  -  (1/X') 

F22  =  (1/YY')  (58) 

F2   =  (1/Y)  -  (1/Y') 

F66  -  (1/S2) 

The  remaining  strength  parameter  F,2  is  obtained  on 
mathematical  grounds.  In  order  for  the  failure  to  be  a  closed 
curve,  the  quadratic  equation  must  be  an  ellipse.  This  leads 
to  the  condition  (see  Ref.(3)  for  details), 

-1  <  F|2/FUF22  <  1 

If  the  Tsai-Wu  criterion  is  required  to  be  a  generalization 
of  the  von  Mises  criterion,  then  (Ft  2/^11^22^  must  be  equal 
to  -0.5.   This  condition  leads  to, 

F12  =  -0.5(FuF22)Va  (59) 

Some  typical  values  of  lamina  strengths  are  presented  in 
Appendix  A. 

10.0  Strength  Ratio  and  Failure  Analysis 

As  discussed  in  the  previous  section,  the  Tsai-Wu  cri- 
terion of  Eq.(56)  can  be  used  to  determine  whether  or  not  a 
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ply  has  failed  by  evaluation  of  the  left  hand  side  of  the 
equation.  By  introducing  a  parameter,  called  the  strength 
ratio,  the  criterion  can  be  used  to  determine  when  a  ply 
fails,  the  sequence  of  ply  failure,  and  the  load  for  laminate 
failure. 

Given  the  temperature  field  of  a  laminate,  T.,  and  the 
relative  magnitudes  of  the  components  of  the  laminate  "load" 
vectors,  N  and  M,  we  seek  to  determine  the  allowable  load 
vectors,  N   and  M  ,  which  will  bring  a  lamina  to  imminent 
failure.  Since  the  laminate  stiffness  matrices,  A  and  D, 
change  with  every  ply  failure,  a  cyclic  (iterative)  pro- 
cedure is  required  to  trace  the  sequence  of  ply  failure. 

The  m-th  cycle  determines  the  m-th  ply  which  will  fail. 
During  the  m-th  cycle,  the  strength  ratio  of  the  i-th  ply  is 
defined  as, 

R.  ,    ,    =  N,  /N,  (60) 

i(m)     la'  1  *■  J 

This  ratio  is  the  same  for  all  components  of  N  and  M,  that 
is,  N.  /N.  =  M.  /M.  ,  for  j=l,2,6.  The  allowable  stresses  in 
a  lamina  during  the  m-th  cycle,  with  respect  to  laminate  axes, 
are  given  by  Eq's.(48), 


(a.  ),    ,=   Q-tA^xCN  ,    .+  N„)  +  z.d;1,(M  ,    ^+  MT) 
^~ia^ (m)    5i  x(m)v~a(m)  ~TJ      i~(m)  ~a(m)   ~T' 


a.6.] 
~i~iJ 


(48) 


where  the  'a1  subscript  denotes  allowable,  the  i  subscript  de- 
notes the  i-th  lamina,  the  m  subscript  denotes  the  m-th  cycle, 

and  A,  -.  and  Dr  -,  are  the  stiffness  matrices  of  the  m-th 
*(m)     ~(m) 

cycle  laminate.  In  Eq. (48) ,  i  takes  on  the  values  of  only  the 
non  failed  lamina.  The  m-th  cycle  laminate  is  comprised  of 
the  lamina  which  have  not  failed  during  the  previous  (m-1) 
cycles . 

It  is  convenient  to  partition  the  lamina  stresses  into 
two  parts;  one  part  due  to  the  N  and  M  loads,  and  another 
part  due  to  the  §T  and  NL.  thermal  vectors.  Denoting  the 
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stresses  due  to  loads  by  a . T ,  and  the  stresses  due  to  temper- 
ature  by  g-T,  we  have, 

a.  r_ ->  =  a .  T  c    \    +  Q  "v  r    ■»  (61) 

-ia(ra)     ~iL(m)     ~iT(m)  *■   ' 

where  the  load  stresses  a.Tr  >,  and  thermal  stresses  a.T(-  ^ 

~iL  (,mj  '  ~il  (.mj 

are  given  by, 

a.,f    ^  =  Q.[  K"fK  -N  r  ^  +   z.Dr\*M  ,  >]     (62) 
~iL(m)     £i L  ~  (m)  ~a(m)      i~(m)  ~a(m)J     *■   v 

a.n  >  =   Q.[  AT1^  "fL  +  z.fi^ii'ftr]  (63) 

~iT(m)     £i   ~  (m)  ~T    i~(m)  ~TJ  v  J 

Using  Eq's.(57),  the  above  lamina  stresses  with  respect  to 
laminate  axes  are  transformed  into  lamina  stresses  with  re- 
spect to  lamina  axes  for  use  in  Eq.(56).  Substituting  lamina 
stresses  with  respect  to  lamina  axes  into  the  Tsai-Wu  cri- 
terion of  Eq.(56)  gives  a  quadratic  equation  for  the  i-th 
lamina  and  m-th  cycle, 

a.r  ^R?r  ^  +  b.,  ^R.r  ->  +  cr  ->  =  0  (64) 

i(m)  i(m)     i(m)  i(m)     l  (m)  v  ' 

where  the  coefficients  are, 


ai(m)=  [FllalL  +  2F12alLa2L  +  F22a2L 

c    2  i 
F66a6LJi(m) 


(65) 


bi(m)=  [2(FllalLalT  +  F12alLa2T  +  F12a2LclT 

+    F22a2La2T    +    F66a6La6T)  (66) 

+    FlalL    +   F2a2lJi(m) 


2  2 

i(m)       L    11    IT  12    IT    2T  22    2T 


+    F66a6T    +    FlalT    +    F2a2T]i(m)    _    1 


24 


(67) 


The  strength  ratios  R.  r    ->  for  each  of  the  lamina  compris- 
ing the  m-th  cycle  laminate  are  obtained  as  the  roots  of  Eq. 
(64) .  The  lamina  with  the  smallest  absolute  strength  ratio  is 
the  lamina  which  will  fail  during  the  m-th  cycle.  Thus  the 
load  factor  to  bring  the  m-th  ply  to  failure  is, 

R*  -v  =  Min[Abs(R.  ,  J]  (68) 

(m)       L    v  i(m)'J  v  * 

Laminate  failure  occurs  when  all  lamina  have  failed.  The 

strength  ratio  (or  load  factor)  for  the  laminate  to  fail,  de- 

noted  by  R ,,,  is  the  maximum  of  the  Rr  >  load  factors,  that  is, 
'   u7  im) 

R  =  Max[R*  .]  (69) 

u      L  (m) J  v  * 

If  the  maximum  strength  ratio  is  obtained  for  m<n,  where  n  is 
the  number  of  plies  comprising  the  laminate,  then  the  lamin- 
ate fails  after  m  plies  have  failed.  The  allowable  loads  that 
the  laminate  can  sustain  without  failure,  for  a  given  temp- 
erature distribution,  are  obtained  from  Eq.(60)  as, 

N   =  R  N 
~a    u~ 


M  =  R  M 
-a    u~ 


(70) 
Illustrative  examples  are  presented  in  section  11. 

11. 0  Results  and  Conclusions 

An  analysis  for  the  prediction  of  the  strength  of  fiber 
reinforced  composite  plates  subjected  to  loads  and  thermal 
environment,  described  in  this  report,  has  resulted  in  a  com- 
puter program  for  implementation  of  the  analysis.  A  listing 
of  the  program,  as  well  as  directions  for  it's  use,  is  given 
in  appendices  B  and  C. 

Several  analyses  have  been  carried  out  in  order  to  com- 
pare the  results  of  the  present  analysis  to  other  analytical 
results  (ref.(4)),  as  well  as  experimental  results.  The  fail- 
ure criterion  used  in  Ref . (4)  is  the  Tsai-Hill  criterion. 
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Figures  6  through  8  present  some  preliminary  results.  Figure 
6  presents  four  sets  of  results  for  cross-ply  laminates  sub- 
jected to  uniaxial  loading: 

1)  experimental  data  for  E-glass/epoxy  from  Ref . (4) , 

2)  analytical  results  for  E-glass/epoxy  from  Ref. (4), 

3)  analytical  results  for  E-glass/epoxy  from  the  present 
analysis,  and 

4)  analytical  results  for  T300/5208  graphite/epoxy  from 
the  present  analysis. 

The  abscissa  M,  called  the  cross-ply  ratio  and  defined  by 

M  =   2tk/  It.    for  even  j,  and  odd  k   (71) 
is  a  measure  of  the  relative  quantity  of  90  and  0  degree 
plies.  It  is  seen  that  for  cross  ply  laminates  under  uniaxial 
loading,  there  is  excellent  agreement  between  experimental 
data  and  analytical  results  for  the  E-glass/epoxy  laminates. 
All  results  are  for  a  curing  temperature  of  270  degrees  Fah- 
renfeit,  and  a  loading  environment  at  70  degrees.  In  each  case, 
the  laminate  did  not  fail  at  first  ply  failure.  The  first  ply 
failure  load  is  given  in  parentheses  below  the  laminate  fail- 
ure load. 

Figure  7  presents  the  same  set  of  results  for  angle  ply 
laminates  as  Figure  6  does  for  cross-ply  laminates.  In  this 
case,  the  laminate  has  the  layup  [-<j>/  +  <|>]  with  equal  thickness 
plies.   It  can  be  seen  that  results  from  the  present  analysis 
for  E  glass/epoxy  laminates  are  in  good  agreement  with  the 
experimental  and  analytical  results  of  Ref.  (4)  for  <j>  less 
than  30°,  and  9  greater  than  60°.   Tsai,  Ref.  (5),  has  a 
forthcoming  report  which  shows  that  the  Tsai-Wu  criterion 
is  in  good  agreement  with  experimental  results  for  T300/5208 
graphite/epoxy  angle  ply  laminates. 

Figure  8  shows  the  effect  of  temperature  on  laminate 
failure  for  t45°  angle  ply  laminates  under  uniaxial  loading 
for  both  E-glass/epoxy,  and  T300/5208  graphite  epoxy  laminates. 

In  order  to  gain  an  understanding  of  the  effect  of 
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temperature  on  laminate  strength,  a  series  of  computer 
analyses  will  be  undertaken.   The  results  of  this  computation- 
al effort  will  be  presented  in  another  report. 
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Appendix  A 

Typical  Properties  for  Two  Lamina 

The  table  below  presents  some  properties  for  two  la- 
mina, E-glass/epoxy ,  and  graphite/epoxy .  The  information 
is  taken  from  references  (3)  and  (4). 


Property 

Units 
SI   (Engr'g.) 

E-gl 

ass/epoxy 

graphite/epoxy 

Ex 

GPa 

(106psi) 

53.8 

(7.8) 

181. 

(26.2) 

E 

y 

V 

xy 
<V 

*   yx 

a 
x 

GPa 

(106psi) 

17.9 

(2.6) 

10.3 

(1.49) 

dimensionless 
GPa   (106psi) 

0.25 
8.6 

(.25) 
(1.25) 

0.28 

7.17 

(.28) 
(1.04) 

dimensionless 

.083 

(.083) 

0.16 

(0.16) 

1/C 

(1/F) 

6.3.y 

(3.5)U 

0.02U 

C.oiDv 

a 

y 

X 

1/C 

(1/F) 

20. S 

(11.4)4i 

22. 5p 

(12.5)u 

MP  a 

(ksi) 

1034. 

(150) 

1500. 

(217.5) 

X' 

MP  a 

(ksi) 

1034. 

(150) 

1500. 

(217.5) 

Y 

MP  a 

Cksi) 

27.6 

(4.) 

40. 

(5.8) 

Y' 

MP  a 

(ksi) 

138. 

(20.) 

246. 

(35.7) 

s 

MP  a 

(ksi) 

41.4 

(6.) 

68. 

(9.9) 

+  vf 

dimensionless 

0.7 

(0.7) 

*  calculated  by  EQ.(12) 

f  volume  fraction  of  fiber 

y  denotes  micro  (i.e., 10   ) 
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Appendix  B 


COMPUTER  LISTING 
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FILE:  STREN1    FORTRAN   A    NAVAL  POSTGRADUATE  SCHOOL 

C       AAA**********  *******************************  A  AAA  A*  A  A  *******  <:  STROO 

C    *                                                            *  STROO 

*     ANALYSIS  OF  THE  STRENGTH  OF  A  SYMMETRIC  COMPOSITE      *  STROO 

C    *       PLATE  DUE  TO  A  THERMAL  ENVIRONMENT                    *  STROO 

C    *                                                              *  STROO 

C       ********  A* ******* ********** A**** AAA*  A  A  AAA AAA  A  AAA  AA  A  A*  A  A  A* A*  *  STROO 

C  STROO 

C    THIS  IS  THE  MAIN  PROGRAM  WHICH  READS  IN  THE  PROPERTIES  OF  THE  STROO 

C    CCPPCSITE,  AS  WELL  AS  THE  INPLANE  LOADS  AND  THERMAL  DISTRIBU-  STROO 

C    TICN  THROUGH  THE  LAMINATE.  IN  ADDITION,  THIS  PROGRAM  CALLS  THE  STROO 

C    SUBROUTINE  CYCLE,  WHICH  IN  TURN  CALLS  SUBROUTINE  FAILUR.  STROO 

C  STROO 

COMMCN    ANGLE (100)  ,Z ( 101)  ,A(3,3 ) ,D(3. 3 ) , AT (3, 3 ) , DT ( 3 ,3 ) , A INV ( 3, 3 ) , WSTROO 

1K( 50),DINV(3.3J ,PLYTHK( 100) ,QBAR11(100) ,  QBAR12U00) ,QBAR22( 100) ,QBSTR0O 

2AR16(100),QBAR26< 100 ) , QB AR66 ( 100), AL PHAl ( 100 ) , ALPHa2( 100) . ALPH A3 ( 1 STROO 

300) ,CELTEM( 100) , STRES1 (100 ) , STRES2 ( 100 ) , STRES3 (100 ) ,TEM( 100 ) , STRESSTROO 

4X(  100), STRSSY(100),TAUXY( 1 00 ),R 1(100) ,R2(100)  , ZN1 , ZN2  ,ZN3 , FXX ,FX ,FSTROO 

5YY,FY,FSS,FXY,QXX,QXY,QYY,QSS,ALFAX,ALFAY,STFRTM,STRSL1(100),  STROO 

6STRSL2C 100) ,STRSL 3(100) , STRSTl(lOO) ,STRST2(100) , STRST3 ( 100 ) ,  STROO 

7STRSTX(100) ,STRSTY( 100) , TAUTXY( 100) ,STRSLX( 100) ,STRSLY(100) ,TAULXYSTROO 

8(1C0),R(70) ,RMAX,ZM1,ZM2,ZM3,  STROO 

9NPLY, ICYCL£,IFAIL»NPL,IT(100)  STROO 

C  STROO 

C         1    AND    2    AFTER    A    NAME    DENOTE    THE    QUANTITY    IS    W.R.T.    THE    PLATE  STROO 

C              AXES    (I.E., THE    GLOBAL    AXES    OF    THE    SYSTEM).  STROO 

C         X    ANC    Y    DENOTE    THE    NATURAL    COORDINATES    OF    A    PLY.    X    IS    ALONG    THE  STROO 

C               FIBER    AND    Y    IS    THE    TRANSVERSE    DIRECTION  STROO 

C         A    ANC    C    ARE    THE    INPLANE    AND    BENDING    STIFFNESS    MATRICES   OF    THE    COMPO-STROO 

C               SITE    PLATE    W.R.T.    LAMINATE    AXES.    THEY    CHANGE    AS    PLIES    FAIL.  STROO 

C         AINV    AND    DINV    ARE    THE     INVERSES    OF    THE    A    AND    D    MATRICES,    RESPECT IVELYSTROO 

C         ALFAX,    AND    ALFAY    ARE   THE   COEFFICIENTS    OF    THERMAL    EXPANSION    W.R.T.  STROO 

C               THE    NATURAL    COORDINATES    OF    A    PLY  STROO 

C         ALPHAKI)  ,ALPHA2(I  )  ,ALPHA3(I  )    ARE    THE    COEFFICIENTS    OF    THERMAL  STRJO 

C              EXPAf^SICN    OF    THE    I-TH    PLY    W.R.T.    THE    LAMINATE    AXES.  STROO 

C         ANGLE(I)     IS    THE    ANGLE    BETWEEN    THE    FIBER    DIRECTION    ( I  .  E.  ,X-AX IS  )    OF  STROO 

C              THE    I-TH    PLY    AND    THE    1-AXIS    OF    THE    LAMINATE    COORDINATE    SYSTEM.  STROO 

C         AT    AND    DT    ARE    THE    A    AND    D    MATRICES   DIVIDED    BY    THE    PLATS    THICKNESS  STROO 

C         DELTEM(I)     IS    THE    TEMPERATURE    DIFFERENCE    OF    THE    I-TH    PLY,    AND  STROO 

C               IS    ECUAL    TO    TEM(  D-STFRTM.  STROO 

C         FXX,FX,FYY,FY,FSS    ARE    THE    FAILURE    PARAMETERS    IN    THE    TSAI-WU    FAILURE  STROO 

C              CRITERION    AND    ARE   OBTAINED    FROM    X,    XPRIME,    Y,    YPRIME,    SS,    AND  STROO 

C              FXYSTR.  STROO 

C         IT(I)    EQUALS    0    FOR    NONFAILED    PLIES,    AND    EQUALS    1    FOR    FAILED    PLIES  STROO 

C         NPLY    IS    THE    NUMBER    OF    PLIES    ABOVE    THE    MIDPLANE    OF    A    SYMMETRIC  STROO 

C               LAMINATE.    THE    TOTAL    NUMBER    OF    PLIES     IS    2*NPLY.    THE    NUMBER    OF  STROO 

C              PLIES    DECREASES    DURING    THE    ANALYSIS    AS    PLIES    FAIL.  STROO 

C         PLYTHK(I)     IS    THE    PLY    THICKNESS    OF    THE     I-TH    PLY.  STROO 

C         QXX,CXY,QYY,AND    QSS    ARE    THE    STIFFNESS    COEFFICIENTS    OF    A    PLY    W.R.T.  STROO 

C              THE    LAMINA    X    AND    Y    AXES.  STROO 

C         QBARll(I)  ,QBAR12(I  ),QBAR22(I  ),QBAR16(I ),QBAR26(I),    AND    QBAR66( I  )  STROO 

C              ARE    THE    STIFFNESS    COEFFICIENTS    OF    THE    I-TH    PLY    W.R.T,    THE    PLATE  STROO 

C               (I.E.,    LAMINATE)     COORDINATES.  STROO 

C         R    IS    THE    LCAD    FACTOR    (I.E., THE    STRENGTH    RATIO)    FOR    THE    GIVEN  STROO 

C              TEMPERATURE    DISTRIBUTION    AND    UNSCALED    INPUT    LOADING.    R    GREATER  STROO 

C              THAN    1.     IS    ASSOCIATED    W.    SAFE    LOADING,    AND    R    LESS    THAN    1.    MEANS  STROO 
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LE:  STREN1 


FORTRAN 


NAVAL  POSTGRADUATE  SCHOOL 


FAILURE  HAS  GCCURRED.  FAILURE 

EQUAL  TO  UNITY. 
STRESKI)  ,STRES2(I  ),STRES3(I  )  AR 

THE  LAMINATE  AXES. 
STPESX(I)  ,STRESY(I)  ,  AND  TAUXYd 

W.R.T.  THE  NATURAL  (PLY)  AXESt 
STRSL  FOLLOWED  BY  A  LAMINA  OR  LA 

STRSL1)  IS  THE  PLY  STRESS  DUE 

COGRCINATES  BEFORE  SCALING  WIT 
STRST  FCLLOWED  BY  A  LAMINA  OR  LA 

DUE  TO  TEMPERATURE  W.R.T.  LAM  I 
STFRTM  IS  THE  STRESS  FREE  TEMPER 
TEM(I)  IS  THE  TEMPERATURE  OF  THE 
THICK  IS  THE  PLATE  THICKNESS 
X,XPRIME,Y,YPRIME,SS,  AND  FXYSTR 

PINEC  EXPERIMENTALLY. 
Z(I)  AND  ZU  +  1)  DENOTE  THE  LOCAT 

SURFACES  OF  THE  I-TH  PLY,  RESP 
ZN1.ZN2,  AND  ZN3  ARE  THE  UNSCALE 
ZNUXY,ZNUYX,EX,EY,  AND  G  ARE  THE 

W.R.T.  PLY  AXES. 


OCCURS  WHEN  ABS ( MIN ( R(  I )  I ) 
E  THE  STRESSES  IN  THE  I-TH 

)  ARE  THE  STRESSES  IN  THE  I 
AFTER  SCALING  WITH  THE  LOA 

MINATE  "SUBSCRIPT"  (I.E.  ST 

TO  LOAD  W.R.T.  LAMINA  OR  LA 

H  THE  LOAD  FACTOR. 

MINATE  "SUBSCRIPT"  IS  THE  P 

NA  OR  LAMINATE  COORDINATES 

ATURE 
I-TH  PLY. 

ARE    THE    FAILURE    PARAMETERS    DETER- 


IQNS    OF    THE    BOTTOM    AND    TOP 
ECTIVELY. 
D    INPLANE    LOADS 
STIFFNESS    COEFFICIENTS    OF 


A    PLY 


PEAC(5,4C)    NPLY,ITYPE,ZNUXY,EX,EY,G, ALFAX,ALF AY, STFRTM 

ZNUYX=ZNUXY*EY/EX 

WRITE    (6,50)    NPLY.ITYPE 


WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 


(6,60) 
(6,70) 
(6,71) 
(6,30) 

(6,90) 
(6,100) 


WRITE(6,105) 

NN=NPLY+1 

NPL^NPLY 


(ZNUXY,ZNUYX, 
ALFAX, ALFAY 
STFRTM 


EX,EY,G) 


IF    ITYPE=0    ALL    PLIES    ARE   THE   SAME    THICKNESS 
IF    ITYPE^l    PLIES   HAVE    DIFFERENT    THICKNESSES 

ICYCLE=1 

Z( 1)=0. 

IF    (  ITYPE.NE.O)    GO   TO    20 

READ    (5,110)    THICK 

PLT=THICK/NPLY 

DO    10    I=1,NPLY 

IT(I)=0 

PLYTHK( I)=PLT 

Z(I+1)=Z(I)+PLT 

CONTINUE 

CONTINUE 

IF    (ITYPE.NE.l)    GO    TO   30 

READ    (5,120)     (J,PLYTHK(I) ,I=1,NPLY) 

DC    25    I=1,NPLY 

Z(I+1)=Z(I)+PLYTHK(I) 

IT(I)=0 

CONTINUE 


IS  STR00560 

STR00570 
PLY  W.R.TSTR00580 
STR00590 
-TH  PLY  STR00600 
D  FACTOR. STR00610 
RSLX  CR  STR00620 
MINATE  STR00630 
STR00640 
LY  STRESSSTR00650 
STR00660 
STR00670 
STR00680 
STR00690 
STR00700 
STR00710 
STR00720 
STRO0730 
STR00740 
STR00750 
STR00760 
STR00770 
STRC0780 
STR00790 
STR00800 
STR00810 
STR00820 
STR00830 
STR00840 
STR00850 
STR00860 
STR00870 
STR00880 
STR00890 
STR00900 
STR00910 
STR00920 
STR009  30 
STR00940 
STR00950 
STR00960 
STR00970 
STR00980 
STR00990 
STR01000 
STR01010 
STR01020 
STR01030 
STR01040 
STRC1050 
STR01060 
STR01070 
STR01080 
STR01090 
STR01100 
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30  READ    (5,130)     ( J , ANGLE ( I) , TEM( I ) , 1=1 , NPLY ) 

WRITE    (6,65) 

WRITE    (6,75)     (J,ANGLE(J),PLYTHK(J) ,TEM( J ),J=1,NPL) 

WRITE    (6,85) 

WRITE    (6,95)     (Z<  I),I=1,NN) 
C 

C  INPLT    THE    APPLIED    INPLANE    LOADS 

C 

READ    (5,140)     ZN1,ZN2,ZN3 

READ(5,140)    ZMi,ZM2,ZM3 

WRITE(6,210) 

WRITE(6,220)    ZN1,ZN2,ZN3 

WRITE(6,230) 

WRITE(6,240)    ZM1,ZM2,ZM3 

WRITE<6,245) 
C 

C         LAMINA    STIFFNESS    COEFFICIENTS    IN    NATURAL    (X,Y)    COORDINATES 
C 

CENCM=1 .0-ZNUXY*ZNUYX 

QXX=EX/DENQM 

QXY=(ZNUXY*EY)/DENOM 

GYY=EY/DENOM 

CSS=G 
C 

WRITE  (6,150) 

WRITE  (6,160)  QXX,«3YY,QXY,QSS 
C 
C    CALCULATE  THE  COEFFICIENTS  OF  THE  FAILURE  CRITERION 


C 


READ    (5,170)    X,XPRIME,Y,YPRIME,SS,FXYSTR 

WRITE(6,175) 

WRITE    (6,180)    X,XPRIME,Y,YPRIME,SS,FXYSTR 

FXX=1./(X*XPRIME) 

FX=(1./X)-(1./XPRIME) 

FYY=1./(Y*YPRIME) 

FY  =  (  l./Y)-< l./YPRIME) 

FXY=FXYSTR*SQRT(FXX*FYY) 

FSS=1./(SS**2) 

WRITE  (6,195) 
WRITE  (6,190) 
WRITE  (6,200)  FXX,FX,FYY,FY,FSS,FXY 

CALL  CYCLE 
RETURN 


C 
C 
40     FORMAT  (2I3,7G10.4l 


50  FORMAT     ( 1H0  ,2X,28HN0.    OF    PLIES    ABOVE    MIDPLANE= , I  3 . 10X,6HI TYPE* , 13) STR01 

60  FORMAT    (  1H0, 2X, 35HI TYPE=0    FOR    UNIFORM    PLY    THICKNESSES)  STR01 

70  FORMAT    ( 3X, 3SHITYPE* 1    FOR    NON    UNIFORM   PLY   THICKNESSES)  STR01 

71  FCRMAT(//1H0,2X, 'LAMINA    ELASTIC    PROPERTIES    W.R.T.    LAMINA    COORDINATSTR01 
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£ES:  •) 
FCRMAT 

1,G9.4,3 
FORMAT 
FCRMAT 
FQRMAT( 
FCRMAT 
FCRMAT 
FORMAT 
FORMAT 
FORMAT 

1) 
FORMAT 

1SS=,G10 
FORMAT 
FORMAT  ( 
FGRMAT 

1RIME=,G 
FCRMAT 

1XY) 
FCPMATi 
FORMAT 
FORMAT( 
FCRMATJ 
FORMAT! 
FORMAT{ 
FOPMAT( 

1Y  THE  S 

2ANC  NCM 
FORMAT 
FORMAT 
FCRMAT 
FORMAT 
END 


( 1H0,2X,5HNUXY=,G9.3, 2X , 5HNUY 

X,2HG=,G9.4) 

(  1H0,2X,7HALPHAX=,G10.4,10X,7 

(//1H0,2X,«THE    STRESS    FREE    TE 

//1H0,2X,' LAMINATE    CONSTRUCTI 

(G10.4) 

(  I3,G10.4) 

( I3.2G10.4) 

(3G10.4) 

(//1H0,2X,46HSTIFFNESS    COEFFI 

(1H0,5X,4HQXX=,G10.3,2X,4HQYY 

.3) 

(6G10.4) 

//1H0,2X.«STRENGTHS    W.R.T.    LA 

( 1H0,4X,2HX=,G10.4,3X,7HXPRIM 

10.4,3X,3HSS=,G10.4,3X,7HFXYS 

{lH0f7X,3HFXX,HX,2HFX,llXt3H 


//lH0t2X 

(4X.G10. 

//1H0,2X 

1H0,5X,» 

//lH0t2X 

1H0,5X,« 

1H0,2X,« 

TRENGTH 

ENTS1  ) 

(1H0,9X, 

( 9X,I3,1 

(1H0.2X, 

(5X.GIO. 


, •TSAI-WU 
4,3X,G10. 
,MNPLANE 
ZN1=» ,G10 
-•MOMENTS 


STRENGTH  PA 
4,3X,G10.4,3 
LOADS  W.R.T 
.4,5X,»ZN2=» 
W.R.T.  LAMI 
ZM1=«  ,G10.4,5X,«ZM2=« 
***  MULTIPL ICATION  OF 
RATIO    AT    LAMINATE    FAI 


STR01660 
X=,G9.3,2X,3HEX=,G9.4,3X,3HEY=STR01670 

STR01680 
HALPHAY=tG10.4)  STR01690 

MPERATURE     IS',G10.4)  STR0L700 

ON:' )  STR01710 

STR0L720 
STR01730 
STR01740 
STR0175Q 
CIENTS    IN    NATURAL    COORD INA TES: STR01760 

STR01770 
*,G10.3,2X,4HQXY=,G10.3,2X,4HQSTR01780 

STR01790 
STR01800 
MINA    AXES: ' )       ,  STR01810 

E=,G10.4,3X,2HY=,G10.4,  3X, 7HYPSTR G1820 
TR=,G10.4)  STR01830 

FYY,  lOXt 2HFY, 10X, 3HFSS, 10X,3HF STRC1840 

STR0185O 
RAMETERS    W.R.T.     LAMINA    AXES:' 
XtG10.4, 3X,G10.4,3X, G10.4) 
.    LAMINATE    AXES:'  ) 
,G10.4,5X, »ZN3=«,G10.4) 
NATE    AXES: * ) 
,G10.4,5X,«ZM3=« ,G10.4) 
THE    ABOVE    LOADS    AND    MOMENTS 
THE    FAILURE    LOADS 


LURE    GIVES 


3HPLY,10X 
OXfG10.4y 
10HZ-LOCA 
4) 


,5HANGLE,11X 
9X,G10.4,6X, 
TION) 


,13HPLY 
G10.4) 


THICKMESSf7X,3HTEM) 


SUBROUTINE    CYCLE 

SUBROUTINE  CYCLE  CONSTRUCTS  THE  A  AND 
FOR  THE  STRESSES  IN  EACH  PLY.  AFTER  TH 
SUBROUTINE  FAILUR  IS  CALLED  TO  OETERMI 
FAVE  FAILED.  SUBROUTINE  FAILUR  IN  TURN 
UNTIL  ALL  PLIES  WHICH  WILL  FAIL  DUE  TO 
ENVIRONMENT    HAVE    FAILED. 


COPMCN 
1K(50),D 
2AR16U0 
300), DEL 
4X(1C0J, 
5YY,FY,F 
6STRSL2( 
7STRSTX( 
8(100), R 
SNPLY,IC 


ANGL 
INV( 
C),Q 
TEM( 
STRE 
SS,F 
100) 
100) 
(70) 
YCLE 


E(100 
3,3)  , 

8AR26 
100)  , 
SY(10 
XY,QX 
,STRS 
,STRS 
,RMAX 
,IFAI 


),Z(101).A(3,3  >,0(3, 
PLYTHK(IOO) ,QBAR11(I 
66( 100),AL 
),STRES2(T 
00) ,R1(100 
QSS,ALFAX, 
RST1  (100), 
UTXY(IOO) , 
M3, 
00) 


( 100) .QBAR 
STRES1Q00 
0),TAUXY(1 
X,QXY,QYY, 
L3(100) ,ST 
TY( 100) ,TA 
,ZM1,ZM2,Z 
L,NPL,IT(1 


)STR01860 
STR01870 
STR01880 
STR01890 
STR01900 
STR01910 
BSTR01920 
STR01930 
STR01940 
STR01950 
STR01960 
STR019  70 
STR01980 
STR01990 
STR02000 
STR02010 
STR02020 
STR02030 
STR02040 
STR02050 
STR02060 
STR02070 
STR02080 
STR02090 
STR02100 

3),AT(3,3) |0T(3f3)tAINV(3t3) ,WSTR 02110 
00) ,  QBAR12U00) ,QBAR22l 100 ) , QBSTR02120 
PHA 1(100) ,ALPHA2(100) , ALPHA3 ( 1STR02130 
00), STRES3(100),TEM( 100) , STRESSTR02140 
),R2(100) ,ZNl,ZN2,ZN3iFXX,FXTFSTR02150 
ALFAY,STFRTM,STRSL1(100) ,  STR02160 

STRST2(100),STRST3<100),  STR02170 

STRSLX(IOO) ,STRSLY(100) ,T AULXYSTR02180 

STR02190 
STR022OO 


D   MATRICES,    AND    THEN    SOLVES 
E    STRESSES    ARE    CALCULATED, 
NE   WHICH,     IF    ANY,    PLIES 
CALLS    SUBROUTINE    CYCLE 
THE    LOAD    AND    THERMAL 
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WRITE    (6,50)     ICYCLE 

NN=NPL+1 
C 
C         ZERO    OCT    THE    A    ANO    D    MATRICESt    AND    THE    THERMAL    LOAD    VECTOR 

A(1,1)=0. 

A(l,2)=0. 
A(  1,3)=0. 
A(2,2Uo. 
A<2,2»»0. 

A(3,3)=0. 

D(1,1)=0. 

D    1,2)=0. 

C(l,3)=0. 

D(2,2)=0. 

C<2,3)=0. 

0<3,3)=0. 

ZNT1=0. 

ZNT2=0. 

ZNT3-0. 

C 

C         LAHINA    STIFFNESS   COEFFICIENTS    W.R.T.    LAMINATE    COORDINATES 

C 

WRITE    (6,100) 
WRITE    (6,110) 

C 

CO    2C    IPLY=1,NPLY 

IF(  IT(IPLY) .EQ.l)    GO    TO    20 

OELTEM(  IPLY)=TEM(  IPLYJ-STFRTM 

ANGLEl=ANGLE(IPLY)/57.2957  8 

S=SIMANGLE1) 

S2=S**2 

S3=S*S2 

S4=S2**2  STR025 

CsCQS(ANGLEl)  STR025 

C2=C**2  STR025 

C3=C2*C  STR025 

C4=C2**2  STR025 

CSSQ=C2*S2  STR026 

(3BAR1K  IPLY)-CXX*C4+2.*(QXY+2„*QSS)*CSSQ+Q\fy*S4  STR026 

QBAR12(  IPLY}=(QXX+QYY-4.*QSS)*CSSQ+QXY*< S4+C41  STR.026 

QBAR22(  IPLY )  =  GXX*S4>2.*(QXY+2.*QSS)*CSSQ+QYY*C4  STR026 

GBAR16( IPLY)=(QXX-QXY-2.*QSS)*S*C3+< QXY-QYY+2. *QSS ) *S3*C  STR026 

CBAR26( IPLY)=(QXX-QXY-2.*QSS )*S3*C  +  (  QXY-QYY+2. *QSS) *S*C3  STR026 

QBAR66( IPLY)=(QXX+QYY~2.*QXY-2.*QSS)*CSSQ+QSS*(S4*C4)  STR026 

C  STR026 

ALPHA1  ( IPLY  )=C2*ALFAX+S2*ALFAY  STR026 

ALPHA2( IPLY)=S2*ALFAX+C2*ALFAY  STR026 

ALPHA3( IPLY|=<ALFAY-ALFAX)*C*S  STR027 

C  STR027 

WRITE    (6,120)     IPLY,QBAR11(  IPLY ) ,QBAR22( I  PLY ) ,QBAR12 ( IPLY ) , QBAR66( I STR 027 

1PLY),0BAR16(IPLY) ,QB AR26< I  PLY )  STR027 

C  STR027 

C    CCNSTRUCT  THE  THERMAL  LOAD  VECTOR,  ZNT  STR027 
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P=PLYTHK(IPLY) 
DTM=CELTEM( IPLY) 
A1=ALPHA1(IPLY) 
A2=ALPHA2(IPLY) 
A2=ALPHA2( IPLY) 


CONSTRUCTICN    OF    LAMINATE    A    AND    D    MATRICES     FROM    LAMINA    MATRICES 

Z1=Z(IPLY+1) 

Z2=Z(IPLY) 

THICK=2.*(Z1-Z2) 

RK*2.*( (Z1**3)-(Z2**3) )/3. 


M1,1)=A(1,1) 
A(l,2)=A(l,2) 
A(  1,2)=A(1.3) 
A(2,1)=A(1,2) 
A(2,2)=A(2,2) 
A(2,2)=A(2,3) 
M3,1)=A(1,3) 
A( 3,2)=A(2,3) 
M3,3)  =  A(3,3)«THICK*QBAR66(IPLY) 


*THICK*QBAR11(IPLY) 
+THICK*CBAR12 (IPLY) 
♦THICK*QBAR16(IPLY) 

♦THICK*Q8AR22 (IPLY) 
♦THICK*QBAR26(IPLY) 


D(  l,l)=D(i,l)+RK*GBARll(IPLY) 

D(lt2)=D(1.2)+RK*CBAR12( IPLY) 

D(l,3)=D(l,3)+RK*C3AR16(iPLY) 

D(2,1)=D(1,2) 

D(2,2)=D(2,2)+RK*GBAR22(IPLY) 

D(2,2)=D(2,3)+RK*QBAR26( IPLY) 

C(3,1)=C(1,3) 

D(3,2)=»D(2,3) 

D(3f 2)=D(3,3)*RK*QBAR66(IPLY) 

CONTINUE 

ZNT1=2.*ZNT1 
ZNT2=2.*ZNT2 
ZNT2=2.*ZNT3 

CALCULATE    THE    TCTAL    INPLANE    FORCES.    TOTNl=  ZNH-ZNT1 

TOTM  =  ZM*ZNTl 
TCTN2=ZN2+ZNT2 
TOTN2=ZN3*ZNT3 


WRITE(6, 125) 
WRITE  (6t130) 
WRITE  (6,140) 
WRITE  (6.150) 

WRITE    (6,155) 


ZN1,ZNT1,T0TN1 
ZN2,ZNT2,T0TN2 
ZN3,ZNT3,T0TN3 


STR02760 
STRC2770 
STR 02780 
STR02790 
STR02800 
STR02810 
STR02820 
TR02830 
TR02840 
STR02850 
STR02860 
STR02870 
STR02880 
STR02890 
STR02900 
STR02910 
STR02920 
STR02930 
STR02940 
STR02950 
STR02960 
STR 02970 
STR029  80 
STR02990 
STR03000 
STR03010 
STR03020 
STR03030 
STR03040 
STR 03050 
STR03060 
STR03070 
STR03080 
STR03090 
STR03100 
STR03110 
STR03120 
STR03130 
STR03140 
ST^03150 
STR03160 
STR02170 
STR03180 
STR03190 
STR03200 
STR03210 
STR03220 
STR03230 
STR03240 
STR03250 
STR03260 
STR03270 
STR03280 
STR03290 
STR03300 
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WRITE  (6,160)  STR033 

DO  25  I=1,NPLY  STR033 

IF(  IT(I).EQ.l)  GO  TO  25  STR033 

WRITE  (6,170)   I,ALPHA1(I) ,ALPHA2( I) ,ALPHA3(I)  STR033 

25     CONTINUE  STR033 

C  STR033 

WRITE  (6,180)  STR033 

WRITE  (6,135)  A(  1,1)  ,A(1,2),A( 1,3) ,D(1,1  ),0(1,2)  ,0(1,31  STR033 

WRITE  (6,190)  A(2,l)  ,A(2, 2), A( 2,3) ,0(2,1) ,0(2,2)  ,0(2,3)  STR033 

WRITE  (6,190)  A(3,l)  ,A(3 ,2 ) , A( 3, 3) ,0 (3, 1 ), D ( 3,  2) ,0(3,3)  STR034 

C  STR034 

C    THE  A  AND  C  MATRICES  PER  UNIT  THICKNESS  OF  LAMINATE  STR034 

C  STR034 

H=2. *Z(NN)  STR034 

WRITE  (6,200)  H  STR034 

AT(1,1)=*A(1,1)/H  STR034 

AT( 1,2)=A( 1,21/H  STR034 

AT(1,3)=A(1,3)/H  STR034 

AT(2,1)*AT( 1,2)  STP034 

AT(2,2)=A(2,2)/H  STR035 

AT(2,3)=A(2,3)/H  STR035 

AT(3tll»ATI 1,3)  STR035 

AT(3,2)*AT(2,3)  STR035 

AT(3,3)=A(3,3)/H  STR035 

C  STR035 

0TU,1)=0(1,1)/H  STR035 

DT(1,2)«C(1,2)/H  STR035 

DT(1,3)»0(1,3)/H  STR035 

DT(2,1)*DT(1,2)  STR035 

CT(2,2)=0(2,2)/H  STR036 

CT(2,3)*D(2,3)/H  STR036 

DT(3,1)=0T(1,3)  STR036 

DT(3,2)=0T(2,3)  STR036 

0T(3,3)=D(3,3)/H  STR036 

C      WRITE  (6,210)  STR036 

C      WRITE  (6,185)  AT ( 1, 1 )  , AT( 1 ,2) , AT( 1 , 3 ) ,DT ( 1 ,1 ) ,DT( 1 , 2) ,0T ( 1 , 3)  STR036 

C      WRITE  (6,190)  AT(2,1  )  , AT (2 ,2 ) , AT( 2, 3  )  ,0T (2, 1 ) , 0T( 2 , 2  )  ,0T( 2, 3 )  STR036 

C      WRITE  (6,190)  AT(3,1) ,AT(3,2) ,AT(3,3| ,DT(3,1) , DT(3 ,2) ,DT(3, 3 )  STR036 


C 


C 


c 


STR036 


C    INVERSION  CF  THE  A  AND  D  MATRICES  STR037 

C  STR037 

CALL  LINV2F  ( A, 3 , 3, A INV, 5 , WK, I ER)  STR037 

CALL  LINV2F  (D,3 ,3,DI NV,5, WK, I ER)  STR037 

C  STR037 

WRITE  (6,220)  STR037 

WRITE  (6,185)  AINV(1,1),AINV(1,2),AINV( 1,3),DINV(1, 1),0INV(  1, 2  ) , DISTR C37 

1NV(1,3)  STR037 

WRITE  (6,190)  AINV(2,1),AINV(2,2),AINV(2,3),0INV(2,1) ,DINV( 2, 2 ) ,DI STR037 

1NV(2,3)  STR037 

WRITE    (6,190)    AINV(3,1) ,AINV(3,2) ,AI NV(3 ,3) ,DI NV( 3 , 1 ) ,DINV ( 3 ,2 ) , DISTR038 

1NV(3,3)  STR038 


STR038 


C         CONSTRUCT    THE    EXPRESSION   FOR    PLY    STRESSES    QUE    TO    TEMPERATURE  STR038 


STR038 


C  1.    AINV*ZNT  STR038 
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WRITE<6,225) 

WRITE(6,230) 

DO  30  I=1,NPLY 

IF(  IT(I  KEQ.l)  GO  TO  30 

AINT1=AINV(  1,1)*ZNT1+AINV(  1, 2) *ZNT2+ AINV {1 ,3 ) *ZNT3 

AINT2=AINV(2,1)*ZNT1+AINV(2,2)*ZNT2*AINV(2,3)*ZNT3 

AINT3*AINV(3,1)*ZNTH-AINV(3,2)*ZNT2*AINV(3,3)*ZNT3 

2.  ALPHA*DELTEM 

ALTM1=ALPHA1(I)*DELTEM<  I) 
ALTM2=ALPHA2U)*DELTEM<  I) 
ALTM2=ALPHA3(I)*DELTEM(I) 

3.  THERMAL    STRESS^    CBAR*{ AI NV*ZNT-ALPHA*DELTEM ) 

A1=AINT1-ALTM1 

A2=AINT2-ALTM2 

A3=AINT3-ALTM3 

STRSTK  I  I  =  QBAR11  (  I)  *A1+QBAR12  (  I  )*A2*QBAR16(  I  )*A3 

STRST2(I)=QBAR12(I)*A1+QBAR22( I  )*A2*QBAR26U  )  *A3 

STRST3(  I1=QBAR16(  I )  *A1+QBAR26(  I  )*A2+QBAR 66< I ) *A3 

WRITE  (6,240) 

CONTINUE 


I ,  STRSTK  I  )  ,STRST2(  I  I  ,STRST3(  I )  ,  DELTEM<  I  ) 


CALCULATE    THE    PLY    STRESSES    DUE    TO    LOADS;    STRSL=Q*AINV*ZN 

DO  35  I=1,NPLY 
IF(ITU).EQ.l)    GO    TO    35 

AINL1=AINV(1,1)*ZNH-AINV(1 ,2 >*ZN2+AI NV( 1, 3 )*ZN3 
AINL2=AINV(2,1)*ZN1+AINV(2,2)*ZN2+AINV(2,3»*ZN3 
AINL3=AINV(3,1)*ZNH-AINV(3,2)*ZN2*AINV(3,3)*ZN3 

STRSLK I)=QBAR11( I) *AINL 1+QBAR 12( I )*A INL2+QBAR16 ( I)*AINL3 
STRSL2C I)=QBAR12(  I) *A  INL1  +  QBAR22( I ) *AINL2+QBAR26< I)*AINL3 
STRSL3<I)=QBAR16(I)*AINLH-QBAR26(I)*AINL2+C8AR66(I)*AIi\IL3 

CCNTINUE 

CALCULATE    THE    STRESSES    DUE    TO    MOMENTS 

IF(Zyi.EC.O..AND.ZM2.EQ.O.  .AND.ZM3 .EQ .0 . )    GO    TO    37 

DO    36    I=1,NPLY 

IF{  IT(I ).EQ.i)    GO    TO    36 

ZZ*<Z<I*l)+Z( I) )/2. 

ZKAP1-DINV<  1,1)*ZM1  +  DINV(1,2)*ZM2*DINVU,3>*ZM3 

ZKAP2=DINV(2,l)*ZMl>CINV(2,2}*ZM2*DINV(2,3i*ZM3 

ZKAP3=DINV<3,1)*ZM1+CINV(3,2)*ZM2+DINV(3,3)*ZM3 

DINMI=ZZ«ZKAP1 
DINM2*ZZ*ZKAP2 
DINM2=ZZ*ZKAP3 


STR03860 
STR03870 
STR03880 
STR03890 
STR03900 
STR03910 
STR03920 
STR03930 
STR03940 
STR03950 
STRC3960 
STR03970 
STR03980 
STR03990 
STR04000 
STR04010 
STR04020 
STR04030 
STR 04040 
STR04050 
STR04060 
STR04070 
STR04080 
STRC4090 
STR04100 
STR04110 
STR04120 
STR04130 
STR04140 
STR04150 
STR04160 
STR04170 
STR04180 
STR04190 
STR 04200 
STR04210 
STR04220 
STR04230 
STR04240 
STRC4250 
STR04260 
STR04270 
STR04280 
STR04290 
STRC4300 
STR 04310 
STR04320 
STR 04330 
STR04340 
STR04350 
STRQ4360 
STR04370 
STR04380 
STR 04390 
STR04400 
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C 

c 
c 

c 


36 
37 


38 
C 
C 
C 

c 


c 

c 
c 
c 


c 

c 
c 


ADO    TC    STRESSES    DUE    TO    INPLANE    LGADS 

STRSL1(II=STRSL1  (D+GBAR11  (I  )*DINM1+QBAR12<  I  )  *DINM2+QBAR  16  (  I) 
S3 

STRSL2( I )=STRSL2( I)*QBAR12( I ) *D INM  H-QBAR 22(  I) *D INM2+QBAR26(  I) 
63 

STRSL3(  I  )=STRSL3(I)+QBAR16(I)*DINM1>QBAR26(I)*DINM2+CBAR66( I) 
£3 

CONTINUE 

CONTINUE 

*RITE(6t270) 

DC    38    I=l,NPLY 

IFUTU ).EQ.l)    GO    TO    38 

WRITE  (6,  280)       I,  STRSLKI)  ,STRSL2(I)  ,STRSL3(I> 

CONTINUE 

FAILURE  ANALYSIS 

9  1.  CONVERT  PLY  STRESSES  TO  NATURAL  (LAMINA)  AXES 

DO  40  I=ltNPLY 

IF(  ITU  ).EQ.l)  GO  TO  40 

ANGLE  WANGLE  ( 1 1/57.29578 

S=SIMANGLE1) 

S2=S**2 

C=C0S(ANGLE1) 

C2=C**2 

SC=S*C 

$TRSTX(I)  =  STRSTliI)*C2-*-STRST2(I)*S2+STRST3m*(2.*SC) 

STRSTYC  I)=STRST1<I)*S2+STRST2(I)*C2-STRST3<I)*(2.*SC) 

TAUTXY( I)=SC*(STRST2(I)-STRST1<I))+STRST3<I)*(C2-S2) 

STRSLXC I)  =  STRSL1( I) *C2+STR SL2( I )*S2+STRS L3(  I )  * (2.* SO 
STRSLY( I )=STRSL1(I)*S2+STRSL2(I)*C2-STRSL3<I)*{2.*SC) 
TAULXYC I)  =  SC*(STRSL2(I)-STRSLl(in*STRSL3<I)*(C2-S2) 

2.  CALCULATE  THE  COEFFICIENTS  OF  THE  QUADRAT IC: A*( R ( I )**2 ) 

♦6*R(I)+C=0. 

AA=FXX*(  STRSLX(I)**2  )+2.*FXY*(  STRSLX  (  I)*STRSLYU)  )+FYY*(STRSL 
1*2)+FSS*(TAULXY( I)**2) 

B=2.*<FXX*STRSLX< I)*STRSTX(I ) +FXY*STRSLX ( I )*STRSTY(I)+ 
1FXY*STRSLY< I ) *STRSTX< I )+FYY*STRSLY ( I )*ST RSTY ( I )+FSS*TAULXY < I) 
2*TAUTXY(  I) )+FX*STRSLX(I)+FY*STRSLY(I  ) 

C=FXX*<  STRSTXU  )**2  J+2.  *FXY*STRSTX  ( I  )*STRSTY  ( I  )  «-FYY*(  STRSTY  (  I 
1*FSS*(TAUTXYI I )**2) +FX*STR STX( I )+FY* STRSTY( I )-l. 

3.  SOLVE    THE    FAILURE    QUADRATIC    FOR    THE    STRENGTH    RATICS*     R( I  ) 

CC=SCRT( ( 8**2 )-4 .*AA*C ) 
Rl(  I)  =  (-B+CC)/(2.*AA) 
R2(I)=(-e-CC)/(2.*AA) 
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STR044 
STR044 
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CONTINUE  STR04960 

STR04970 
WRITE    <6,250)  STR04980 

00    4*-i~"t,NPLY  ••  STR04990 

IF(  IT(I  ).EQ.l)    GO    TO    41  STR05000 

*RITE    (6,260)       IiSTRSTXCI)  ,STRSTYU)  ,TAUTXY(II  ,Rim,R2(II  STR05010 

CONTINUE  STR05020 

STR05030 
WRITE(6,290)  STR05040 

CC    42    I=1,NPLY  STR05050 

IF(  IT(I J.EQ.l)    GO    TO    42  STR05060 

WRITE(6,300)       ItSTRSLX(I) ,STRSLYUI ,TAULXYUI  STR05070 

CONTINUE  STR0508O 

STR05090 
CALL    FAILUR  STR05100 

STR05110 
RETURN  STR05120 

STR05130 
FORMAT    <1H1,1X,7HICYCLE*,I3)  STR05140 

FORMAT  (//1H0,2X,«PLY  STIFFNESS  COEFFICIENTS  W.R.T.  LAMINATE  AXES: STR05 150 
1«)  STR05160 

FORMAT  (1H0,3X,3HPLY,3X,6HQBAR11,5X,6HQBAR22,5X,6HQBAR12,5X,6HQBARSTR05170 
166,5X,6HQBAR16,5X,6HQBAR26I  STR05180 

FORMAT    <2X,I3,6<  1X,G10.3) )  STR05190 

F0RMAT(//1H0,2X, 'APPLIED,    THERMAL*    AND   TOTAL    INPLANE    FORCES:1)  STR05200 

FORMAT  (1H0,5X,4HZN1=,G10. 4, 10X, 5HZNT1= ,G10.4 , 10X ,6HTCTN1= ,G10.4)  STR05210 
FCRMAT     (6X,4HZN2=,G10.4,10X,5HZNT2=,G10.4,10X,6HT0TN2=,G10.4)  STR05220 

FORMAT    ( 6X,4HZN3=,G10.4,10X,5HZNT3=,G10.4,10X,6HTOTN3=,G10.4)  STR05230 

FORMAT  (//1H0.2X, ^THERMAL  COEFFICIENTS  OF  EXPANSION  W.R.T.  LAMINATSTR05240 
IE    AXES:* )  STR05250 

FORMAT    ( 1H0,3X,3HPLY,11X, 6HALPHA1 , 14X.6H ALPHA2 , 14X  ,6FALPHA6 )  STR05260 

FORMAT    <3X,I3,10X,G10.4,10X,G10.4, 10X.G10.4)  STR05270 

FORMAT    (//1H0,12X,8HA-MATRIX,35X,8HD-MATRIX)  STR0528O 

FORMAT  (1H0,  IX,  Gil.  4,  IX,  Gil.  4, 1X,G11 .4, 8X,311  .4, IX ,G11. 4, IX ,G11 .4) STR05290 
FCRMAT    (2X,G11.4,lXtG11.4,lX.G11.4,8X,G11.4, IX »G1 1. 4, 1X,G11 . 4)  STR 05300 

FORMAT    (//1H0,5X,16HPLATE    THICKNESS* ,G10 .4)  STR05310 

FCRMAT  </lH0,9X,18HA-MATRIX/THICKNESS,25X,18HD-MATRIX/THICKNESS)  STR053  20 
FORMAT    (1H0,9X,16HA-INVERSE    MATRIX, 30X.16HD-INVERSE    MATRIX)  STR05330 

FORMAT( //1H0,2X,»STRESSES  IN  LAMINA  AND  LAMINATE  COORDINATES  BEF0RSTR05340 
£E    MULTIPLICATION    WITH   THE    LOAD    FACTOR:*)  STR05350 

FORMAT  (1H0,3X,3HPLY,6X,6HSTRST1,12X,6HSTRST2,14X,6HSTRST3,9X,6HDESTR 05360 
1LTEM)  STR05370 

FCRMAT    (4X,I3,7X,G10.4,8X,G10.4,9X,G10.4,5X,G10.4)  STR05380 

FORMAT  ( 1H0,3X,3HPLY,6X,6HSTRSTX,9X,6HSTRSTY,1IX,6HTAUTXY, 12X , 2HR1STR05390 
1,12X,2HR2)  STR05400 

FORMAT    (4X,I3,5X,G10.4,5X,G10.4,5X,G10.4,5X.G10.4,5X,G10.4)  STR 054 10 

FORMAT! 1H0,3X,3HPLY,6X,6HSTRSL1,12X. 6HSTRSL2 ,14X,6HSTRSL3 I  STR05420 

FCRMAT(4X,I3,7X,G10.4,8X,G10.4,9X,G10.4)  STR 05430 

FORMAT! 1H0,3X,3HPLY,6X,6HSTRSLX,9X26HSTRSLY,11X,6HTAULXY)  STR 05440 

F0RMAT(4X,I3,5X,G10.4,5X,G10.4,5X,G10.4)  STR05450 

END  STR05460 

STR05470 

STR05480 
SUBROUTINE    FAILUR  STR05490 

STR05500 
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C    THIS  SUBROUTINE  DETERMINES  WHICH  PLIES  HAVE  FAILED  BY  USING 

C    THE  TSAI-WU  FAILURE  CRITERION.  THE  PLY  WHICH  HAS  FAILED  IS 

C    ELIMINATED  AND  THE  ANALYSIS  IS  RETURNED  TO  SUBROUTINE  CYCLE, 

C    WHICH  RECALCULATES  THE  STRESSES  FOR  THE  MODIFIED  LAMINATE  AND 

C    RETURNS  TO  SUBROUTINE  FAILUR  AGAIN.  THIS  PROCESS  IS  REPEATED 

C    UNTIL  ALL  PLIES  WHICH  ARE  TO  FAIL  HAVE  FAILED. 

C 

C 

C         THIS    SLBROUTINE    CHECKS    FOR    PLY    FAILURE. 

C 

COMMCN  ANGLE ( 100) .Z(101)  ,A(3,3)  ,D( 3 , 3  ) , AT( 3 ,3)  ,DT ( 3 . 3 ) , AINV ( 3, 3 
lK(50)fDINV(3,3),PLYTHK(100),QBARll( 100) ,QBAR12(100) ,CBAR22(100) ,QBSTR056 
2AR16(1QO),QBAR26(100),QBAR66(100),ALPHA1 (100 ) , ALPHA2( 100) , ALPHA3( 1STR05* 
300)  ,DELTEM(  100) , STRES1 ( 100 ) , STRES2( 100)  .STRES3  (1001  ,TEM(  1.00  ) .  STRESSTR054 
4X(100),STRESY(100),TAUXY( 1 00 ),R 1(100) ,R2(100) ,ZN1 ,ZN2 ,ZN3 ,F XX ,FX ,FSTR05* 
5YY, FY, FSS, F XY, QXX  , QXY, QYY , QS S , ALF AX, ALFA Y, STFRTM , STRSL 1(100),  STR05* 

6STRSL2( 100) .STRSL3U00)  ,  STRST1  ( 100) , STRST2(100) ,STRST3(100) ,  STR05* 

7STRSTX(100),STRSTY(100),TAUTXY(100) , STRSLX( 100 ), STRSL Y( 100) , TAULXYSTR05* 
8(100),R(70) ,RMAX,ZM1  ,ZM2,ZM3, 
SNPLY,ICYCLE,IFAIL,NPL,IT( 100) 
C 

R(  ICYCLE  Ml. E+20 


10 

c 

c 
c 
c 

DO    10    I=1,NPLY 

IF(  IT(I J.EQ.l)    GO    TO    10 

A8SR1=A£S(R1(I) ) 

ABSR2=ABS(R2(I) ) 

ABSR=AMIN1(ABSR1,ABSR2) 

IF(ABSR.LT.R( ICYCLE) )    IFAIL=I 

IF (ABSR.LT.Rt ICYCLE)  )    R< ICYCLE )=ABSR 

CONTINUE 

WRITE i6, 110)     ICYCLE, I FAIL,R( ICYCLE) 

CALCULATE    THE    STRESSES    AT    IMPENDING    PLY    FAILURE 

CC    20    I=1,NPLY 

IF( ITC I ).EQ.l)    GO    TO    20 

STRESX( I)=R(ICYCLE)*STRSLX(II+STRSTX(I) 

STRESY( I)=RCICYCLE)*STRSLY (I )+STRSTY( I) 

TAUXVCI )=R(ICYCLE)*TAULXY( I)*TAUTXY(  I) 
20  CONTINUE 

C 

WRITE(6,100) 

WRITE(6, 120) 

DC  3C  I=1,NPLY 

IF(  ITC I l.EQ.l)  GO  TO  30 

WRITE (6, 125)  (I,STRESX(I),STRESY(I) ,TAUXY(I) ) 
30     CONTINUE 
C 
C    CALCULATE  THE  STRENGTH  RATIO  FOR  LAMINATE  FAILURE 


C 


IFUCYCLE.EC.il    RMAX=R(1) 
RM AX  =AMAX1(R(  ICYCLE  ),RMAX) 
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CHECK  FOR  LAMINATE  FAILURE  STR06060 

STR06070 
NPl=NPL-l  STR06080 

IF(NPL.NE.O)  GO  TO  50  STR06090 

WRITE(6,130)  RMAX  STR06100 

IF(NPL.EC.O)  STOP  1  STR06110 

ICYCLE  =ICYCLE-H  STR06120 

IT(IFAIL)=1  STR06130 

STR06140 
CALL  CYCLE  STR06150 

STR06160 
RETURN  STP06170 

STR06180 

0         F0RMAT(//1H0,2X,« STRESSES    AT     IMPENDING    PLY    FAILURE    AFTER    MULT  IPL ICSTR06190 

SATION    W.    THE    LOAD    FACTOR:')  STR06200 

0         FORMAT (//lH0f2X, • CYCLE=« , I 3, 5X , • IFAI L=» , 13, 5 X, 'LOAD    FACTOR*1 ,G 10.4STR06210 

1)  STR06220 

0         FCRMAT(  1H0,3X, «PLV«,8X, • STRESX' ,8X,  •  STRESY* ,  10X, • JAUXY' )  STR06230 

5         FGRMAT(4X,I3,6X,G10.4,4Xf G10.4,4X,G10.4)  STR06240 

0         FORMAT< ///////1H0,2X,»***       THE    STRENGTH    RATIO    FOR    THE    LAMINATE    IS«STR06250 

£fG10.4)  STR06260 

STR06270 
END  STR06280 
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Appendix  C 

Instructions  for  Use  of  the  Computer  Program 

The  following  describes  the  sequence  of  input  data  for  the 
computer  program: 

1.  Read  (5,40)  NPLY,  ITYPE,  ZNUXY,  EX,  EY,  ALFAX,  ALFAY, 

STFRTM 
40  Format  (213,  7G10.4) 

NPLY  =  Number  of  plies  above  the  laminate  midplane 

ITYPE  is  0  if  all  plies  are  of  equal  thickness; 

ITYPE  is  1  if  plies  have  different  thicknesses 

ZNUXY  (v   )  =  lamina  longitudinal  Poissons  ratio 
xy 

EX  (Ex)  =  lamina  longitudinal  stiffness 

EY  (Ey)  =  lamina  transverse  stiffness 

ALFAX  (a  )  =  lamina  longitudinal  coefficient  of 

thermal  expansion 
ALFAY  (a  )  =  lamina  transverse  coefficient  of  thermal 

yJ 

expansion 
STFRTM  =  stress  free  temperature  (i.e.,  the  curing 
temperature  of  the  laminate) 

2.  Read  (5,110)  THICK  (not  inputted  if  ITYPE  =  1) 
110  Format  (G10.4) 

THICK  =  laminate  thickness 

3.  Read  (5,120)  (J,PLYTHK  (I),  I  =  1,  NPLY)  (not  inputted 

if  ITYPE  =  0) 
120  FORMAT  (13,  G10.4) 
J  =  ply  number 
PLYTHK  (I)  =  ply  thickness 

4.  Read  (5,130)  (J,  ANGLE  (I),  TEM  (I),  I  =  i,  NPLY) 
130  Format  (13,  2G10.4) 

J  =  ply  number 

ANGLE  (I)  =  <$>  •    =  ply  orientation  angle 

TEM  (I)    =  T.  =  ply  temperature 


46 


5.  Read  (5,140)  ZN1,  ZN2,  ZN3 
140  Format  (3610.4) 

ZN1  =  N,  =  normal  force  along  laminate  1-axis 
ZN2  -  N2  =  normal  force  along  laminate  2-axis 
ZN3  =  N-  «  shear  force  on  laminate 

6.  Read  (6,140)  ZM1 ,  ZM2 ,  ZM3 

ZM1  =  M,  =  bending  moment  on  laminate  1-face 
ZM2  =  M-  =  bending  moment  on  laminate  2- face 
ZM3  =  M,  =  twisting  moment  on  laminate 

7.  Read  (5,170)  X,  XPRIME,  Y,  YPRIME,  SS,  FXYSTR 
170  Format  (6G10.4) 

X  =  longitudinal  tensile  strength  of  lamina 

XPRIME  (X')  =  longitudinal  compressive  strength  of  lamina 

Y  =  transverse  tensile  strength  of  lamina 

YPRIME  =  (Y1)  =  transverse  compressive  strength  of  lamina 

SS  (S)  =  lamina  shear  strength 

FXYSTR  (F*  )  =  O.S  (in  accordance  with  Eq.  (59)) 
xy 

A  sample  output  of  a  typical  analysis  is  presented  on  the 
pages  following  this  page.   This  is  the  problem  of  a  cross-ply 
E-glass/epoxy  laminate  subjected  to  uniaxial  loading,  i.e. 
(N,,  N2,  N-)  =  (1,0,0).   The  cross-ply  ratio,  M,  is  0.2.   The 
curing  temperature  is  270  F,  and  the  temperature  of  the  lamin- 
ate is  70°F. 
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